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ABSTRACT 

The  work  accomplished  during  the  second  year  of  this  research  project  is 
combination  of  analytical  and  experimental  investigations: 


t 

a 


The  analytical  part  deals  with  the  nonlinear  response  of  a  three-degree-of- 
freedom  aeroelastic  structural  model  in  the  neighborhood  of  combination 
internal  resonance  condition.  The  Fokker-Planck  equation  approach  is  used 
to  derive  a  general  differential  equation  for  the  response  statistical  joint 
moments.  This  equation  is  found  to  constitute  a  set  of  infinite  coupled 
first  order  differential  equations.  In  view  of  the  system  complexity  an 
attempt  is  made  to  close  the  infinite  hierarchy  by  using  a  Gaussian  scheme. 
This  scheme  leads  to  27  differential  equations  in  the  first  and  second 
response  moments.  The  equations  are  solved  by  using  numerical  integration. 
The  solution  shows  that  the  response  coordinates  are  non- stationary  random 
processes  and  the  three  normal  modes  are  in  complete  nonlinear  interaction. 
The  interaction  is  found  to  be  very  strong  at  a  region  of  internal  detuning 
which  is  shifted  from  the  exact  internal  resonance  condition.  This  result 
is  under  further  investigation  by  using  a  non-Gaussian  closure  scheme. 


The  experimental  investigation  is  conducted  out  on  a  two  degree -of- freedom 
model  whose  analytical  solution  was  obtained  during  the  first  year  of  this 
project.  When  the  first  normal  mode  is  externally  excited  by  a  band  limited 
random  excitation,  the  system  mean  square  response  is  found  to  be  linearly 
proportional  to  the  excitation  spectral  density  level  up  to  a  certain  level 
above  which  the  two  normal  modes  exhibit  discontinuity  governed  mainly  by 
the  internal  detuning  parameter  and  the  system  damping  ratios.  The  results 
are  completely  different  when  the  second  normal  mode  is  excited.  For  small 


levels  of  excitation  spectral  density  the  response  is  dominated  by  the 
second  normal  mode.  For  higher  levels  of  excitation  spectral  density  the 
first  normal  mode  attends  and  interacts  nonlinearly  with  the  second  mode  in 
a  form  of  energy  exchange . 

New  directions  of  this  research  project  have  been  evolved  during  this  year. 
These  include  the  influence  of  random  in  aerodynamic  forces  on  the  nonlinear 
response  of  typical  aeroelastic  structures.  Two  aeroelastic  models,  are 
chosen  to  carry  out  this  investigation 


INTRODUCTION 


This  report  presents  the  main  results  of  the  research  project  "Stochastic 
Nonlinear  Flutter  of  Aeroelastic  Structures"  funded  by  the  AFOSR  under  grant 
No.  AFOSR- 85- 0008 .  The  report  covers  only  the  work  performed  during  the 
second  year  of  this  project.  Furthermore,  additional  new  research  problems 
have  been  evolved  during  this  year.  These  problems  include  the  influence  of 
the  aerodynamic  forces  on  the  random  response  of  nonlinear  aeroelastic 
structures.  A  formal  proposal  will  be  submitted  next  June  for  new  funds  to 
support  these  new  problems. 

SUMMARY  OF  MAIN  RESULTS 

ANALYTICAL  INVESTIGATION 

The  linear  and  autoparametric  modal  interactions  in  a  three  degree-of- 
freedom  structure  subjected  to  wide  band  random  excitation  are  examined. 
For  a  structure  with  constant  parameter  properties  the  linear  response  is 
obtained  in  a  closed  form.  When  the  structure  stiffness  matrix  involves 
random  fluctuations,  the  governing  equations  of  motion,  in  terms  of  normal 
coordinates,  are  found  to  be  coupled  through  parametric  terms.  The  struc¬ 
ture  response  is  mainly  governed  by  the  condition  of  mean  square  stability. 
The  boundary  of  stable -unstable  response  is  obtained  as  a  function  of  the 
internal  detuning  parameter.  The  results  of  the  linear  system  with  constant 
parameters  are  used  as  a  reference  to  measure  the  deviation  of  the  system 
response  when  the  nonlinear  inertia  coupling  is  included.  In  the  neighbor¬ 
hood  of  combination  internal  resonance  the  system  random  response  is  deter¬ 
mined  by  using  the  Fokker- Planck  equation  approach  together  with  the 
Gaussian  closure  scheme.  This  approach  results  in  27  coupled  first  order 
differential  equations  in  the  first  and  second  response  moments.  These 


equations  are  solved  by  numerical  integration.  The  response  is  found  to 
deviate  significantly  from  the  linear  solution  when  the  system  internal 
detuning  is  close  to  the  exact  internal  resonance.  The  autoparametric 
interaction  is  found  to  depend  significantly  on  the  system  damping  ratios 
and  the  nonlinear  coupling  parameter.  In  the  vicinity  of  combination  inter¬ 
nal  resonance,  the  second  normal  mode  mean  square  exhibits  an  increase 
associated  with  a  corresponding  decrease  in  the  first  and  third  normal 
modes.  The  first  normal  mode  shows  a  very  small  deviation  from  the  linear 
solution  which  implies  that  the  nonlinear  interaction  takes  place  between 
the  second  and  third  normal  modes .  This  unexpected  feature  is  currently 
under  further  investigation  in  parallel  to  a  non-Gaussian  closure  analysis. 
The  results  of  this  work  has  been  accepted  for  publication  in  the  journal  of 
Probabilistic  Engineering  Mechanics.  A  copy  of  the  page  proof  of  this  paper 
is  attached. 

In  order  to  enhance  our  understanding  to  the  main  results  of  structural 
dynamics  with  parameter  uncertainties,  the  P.I.  has  conducted  an  extensive 
literature  survey  which  has  been  accepted  for  publication  in  che  ASME 
Applied  Mechanics  Reviews,  A  preprint  of  this  paper  is  attached. 

EXPERIMENTAL  INVESTIGATION 

A  series  of  experimental  tests  is  conducted  on  a  two  degree-of- freedom 
elastic  structural  model.  The  model  is  subjected  to  a  band-limited  random 
excitation  with  a  central  frequency  very  close  to  one  of  two  normal  mode 
frequencies.  The  band  width  is  selected  such  that  only  the  mode  under 
consideration  is  excited.  The  model  normal  mode  frequencies  are  adjusted  to 
have  the  ratio  2  to  1 .  This  ratio  meets  the  condition  of  internal  resonance 
of  the  analytical  model.  When  the  first  normal  mode  is  external  excited  the 


system  mean  square  response  is  found  to  be  linearly  proportional  to  the 
excitation  spectral  density  up  to  a  certain  level  above  which  the  two  normal 
modes  exhibit  discontinuity  governed  mainly  by  the  internal  detuning 
parameter  and  the  system  damping  ratio.  The  results  are  completely 
different  when  the  second  normal  mode  is  externally  excited.  For  small 
levels  of  excitation  spectral  density  the  response  is  dominated  by  the 
second  normal  mode.  For  higher  levels  of  excitation  spectral  density  the 
first  normal  mode  attends  and  interacts  with  the  second  normal  mode  in  a 
form  of  energy  exchange.  A  number  of  deviations  from  theoretical  results 
are  observed  and  discussed  in  the  attached  manuscript  (28th  SDM  Conference, 
Paper  No.  87-0079-CP)  which  will  be  presented  at  the  AIAA/ASME/ASCE/AHS  28 
Structures,  Structural  Dynamics  and  Materials  Conference. 

NEW  RESEARCH  DIRECTIONS 

The  influence  of  random  aerodynamic  forces,  in  subsonic  and  supersonic  flow 
regimes,  on  the  nonlinear  response  of  typical  aroelastic  structural  models 
has  been  identified  as  a  potential  problem  in  aeroelastic  flutter.  Two 
models,  which  include  a  cantilovel  wing  beam  and  a  flat  panel,  will  be 
adopted.  The  equations  of  motion  of  the  first  model  have  been  derived  by 
using  the  Lagrangian  formulation.  The  aerodynamic  forces  are  derived  by 
using  the  modified  strip  theory  which  includes  the  effect  of  the  span  finite 
length.  The  linear  part  of  the  equations  of  motion  has  ben  considered  to 
derive  the  flutter  boundaries  and  to  identify  the  system  eigenvalues.  This 
prelimenary  analysis  is  essential  to  identify  the  system  parameters  which 
satify  the  condition  of  internal  resonance  between  bending  and  torsional 
motions.  The  work  of  this  part  is  in  progress. 
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Stochastic  modal  interaction  in  linear  and  nonlinear 
aeroelastic  structures 


R.  A.  Ibrahim  and  Z.  Hcdayati* 

Texas  Tech.  University.  Department  of  Mechanical  Enifinccnnif,  Luhhoik.  Tc\us  79409.  USA 

The  linear  and  auloparamctric  modal  interactions  in  a  three  defree-of-freedom  structure  under 
wide  band  random  excitation  are  examined.  For  a  structure  with  constant  parameters  the  linear 
response  is  obtained  in  a  closed  form.  When  the  structure  stiffness  matrix  involves  random 
fluctuations,  the  governing  equations  of  motion,  in  terms  of  the  normal  coordinates,  are  found  to 
be  coupled  through  parametric  terms.  The  structural  response  is  mainly  governed  by  the  condition 
of  mean  square  stability.  The  boundary  of  stable-unstable  responses  is  obtained  as  a  function  of 
the  internal  detuning  parameter.  The  results  of  the  linear  system  with  constant  parameters  are 
used  as  a  reference  to  measure  the  deviation  of  the  system  response  when  the  nonlinear  inertia 
coupling  is  included.  In  the  neighbourhood  of  combination  internal  resonance  the  system  random 
response  is  determined  by  using  the  Fokker  Planck  equation  approach  together  with  the  Gaussian 
closure  scheme.  This  approach  results  in  27  coupled  first  order  differential  equations  in  the  first 
and  second  response  moments.  These  equations  are  solved  numerically.  The  response  is  found  to 
deviate  significantly  from  the  linearsolution  when  the  system  internal  detuning  is  close  to  the  exaci 
internal  resonance.  The  autoparametric  interaction  is  found  to  depend  significantly  on  the  system 
damping  ratios  and  a  nonlinear  coupling  parameter.  In  the  vicinity  of  combination  internal 
resonance,  the  second  normal  mode  mean  square  exhibits  an  increase  associated  with  a 
corresponding  decrease  in  the  first  and  third  normal  modes. 


J.  INTRODUCTION 

The  modal  analysis  of  acrolastic  structures  is  usually 
carried  out  by  using  one  of  the  available  computer  codes 
for  eigenvalues  and  eigenvectors.  These  computer 
algorithms  are  useful  in  determining  the  structural 
dynamic  behavior  under  various  types  of  excitations.  The 
fust  step  usually  involves  the  detei  initiation  of 
eigenvalues  and  eigenvectors.  With  this  information  one 
can  determine  the  linear  response  to  deterministric  or 
random  excitations.  For  systems  with  constant 
parameters  the  mean  square  response  to  external  white 
noise  is  linearly  proportional  to  the  excitation  spectral 
density.  If  the  excitation  is  acting  parametrically  to  the 
system  the  equilibrium  state  could  be  stable  or  unstable  in 
a  stochastic  sense.  In  certain  situations  the  structure  may 
not  behave  according  to  the  linear  theory  of  small 
oscillations  and  a  number  of  complex  response 
characteristics  such  as  amplitude  jump,  internal 
resonance,  saturation  phenomenon,  and  chaotic 
motion1,2  may  be  observed  These  new  characteristics 
owe  their  origin  to  the  system  inherent  nonlinearitics 
which  should  not  be  ignored  in  dynamic  analysis 
In  aircraft  structures  several  types  of  nonlinearitics 
have  been  reported.  Brcitbach3  classified  structural 
nonlinearitics  into  distributed  and  concentrated 
Distribution  non'incarity  is  induced  by  elastic 
deformation  in  riveted,  screwed  and  bolted  connections 

•  Current!)  PhD  Muiicnl.  M.tss;ichiKCU>  Insiikiie  of  1  ccImw'Iojv 
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as  well  as  within  the  structural  components  themselves. 
Concentrated  nonlinearity  acts  locally  lumped  in  control 
mechanisms  or  in  the  connecting  parts  between  wing  and 
external  stores.  This  nonlinearity  results  from  back-lash 
in  the  linkage  elements  of  the  control  system,  dry  friction 
in  control  cable  and  push  rod  ducts,  kinematic  limitation 
of  the  control  surface  deflection,  and  application  of  sprine 
tab  system  provided  for  relieving  pilot  operation. 
Breitbach4  determined  the  flutter  boundaries  for  three 
different  configurations  distinguished  by  different  types  of 
nonlinearitics  in  the  rudder  and  aileron  control  system  of 
a  sailplane  It  was  shown  that  the  influence  of  hvsteretic 
damping  results  in  a  considerable  stabilizing  effect  and  an 
increase  in  the  flutter  speed .  However,  this  special  type  of 
non-linearity  does  not  bring  the  structural  response  into  a 
bounded  limit  cycle.  Similar  elTects  of  nonlinearitics  due 
to  friciion  and  back -lash  w  ere  considered  by  De  Ferrari  et 
ul Peloubet  ci  al ",  Reed  et  al.  and  Desmarais  and 
Reed8  examined  the  effects  of  control  system 
nonlinearitics,  such  as  actuator  force  or  deflection  limits, 
on  the  performance  of  an  active  flutter  suppression 
system.  It  was  shown  that  a  nonlinear  system  which  is 
stable  with  respect  to  small  disturbances  may  be  unstable 
with  respect  10  large  ones  Another  important  feature  was 
that  a  stoic  on  a  pylon  with  low  pitch  stiffness  can  provide 
substantial  increase  m  flutter  speed  and  reduce  the 
dependency  of  flutter  on  the  mass  and  inertia  of  stores 
relative  to  that  of  stiff-mounted  stores. 

In  stiuctural  dynamics,  the  nonlinearity  may  take  one 
of  three  classes'*  l0:  elastic,  inertia  and  damping 
nonlinearitics.  Hlaxtic  nonlinearity  stems  from  nonhneai 


strain-displacement  relations  which  are  inevitable.  Inertia 
nonlinearity  is  derived,  in  Lagrangian  formulation,  from 
the  kinetic  energy.  In  multi-degree-of-freedom  systems 
the  normal  modes  may  involve  nonlinear  inertia  coupling 
which  may  give  rise  to  what  are  effectively  parametric 
instability  phenomena  within  the  system.  The  parametric 
action  is  not  due  to  the  external  loading,  as  in  the  case  of 
parametric  vibration,  but  to  the  motion  of  the  system 
itself  and,  hence,  is  described  as  autoparametric11.  The 
main  feature  of  autoparametric  coupling  is  that  responses 
of  one  component  of  the  structure  give  rise  to  loading  of 
another  component  through  time-independent 
coefficients  in  the  corresponding  equation  of  motion.  The 
deterministic  autoparametric  interactions  in  two  and 
three  freedom  systems  were  examined  by  Barr  and 
Ashworth12,  Haddow  et  a! .'3,  Ibrahim  et  al .'4,  and 
Ibrahim  and  Woodal15.  These  studies  have  shown  that 
the  mode  which  is  externally  excited  exhibits  a  saturation 
phenomenon  in  which  energy  is  transferred  to  other 
modes  involved  in  the  nonlinear  coupling.  The  stochastic 
aspects  of  parametric  and  autoparametric  vibrations  have 
recently  been  documented  in  a  recent  research 
monograph  by  Ibrahim16. 

To  the  authors’  knowledge  the  random  response  of 
systems  with  autoparametric  coupling  has  been  restricted 
to  two-degree-of-freedom  systems.  This  paper  deals  with 
the  linear  and  nonlinear  modal  interactions  of  a  three 
degree-of-freedom  aeroelastic  structure  subjected  to 
random  excitation.  The  deterministic  respinses  of  this 
model  under  various  internal  resonance  conditions 
Y,k,<j;=  0  (where  kt  are  integers  and  cu,  are  tge  system 
normal  mode  frequencies)  have  been  determined  by 
Ibrahim  et  a/.14,15.  The  system  involves  quadratic 
nonlinear  inertia  which  couples  the  system  normal 
modes.  It  was  shown  that  under  principal  internal 
resonance,  the  mode  which  is  directly  excited  is 
suppressed  and  energy  is  transferred  to  the  other  mode. 
When  the  structure  possess  combination  internal 
resonance  of  the  summed  type  the  normal  mode 
amplitudes  did  not  achieve  a  steady  stale  and  the 
response  is  characterized  by  energy  exchange  between  the 
three  modes. 

The  main  objectives  of  this  paper  are  to  present  the 
linear,  parametric  and  autoparametric  random  responses 
of  the  same  aeroelastic  model  considered  in  Refs  14  and 
15.  The  mean  square  responses  will  be  evaluated  for  a 
model  with  constant  parameters  and  for  a  model  with 
random  variations  in  its  stiffness  matrix.  The  nonlinear 
random  response  of  the  system  in  the  neighbourhood  of 
combination  internal  resonance  of  the  summed  type  will 
be  determined  by  using  the  Fokker  Planck  equation 
approach  together  with  a  Gaussian  closure  scheme.  The 
effects  of  the  system  nonlinearity  and  damping 
coefficients  on  the  mean  square  responses  will  be 
examined. 


II.  BASIC  MODEL  AND  EQUATIONS  OF 
MOTION 

Fig.  1  shows  a  schematic  diagram  of  an  analytical  model 
of  an  aircraft  subjected  to  random  excitation  F(t).  The 
fuselage  is  represented  by  the  main  mass  m3,  linear  spring 
K j,  and  dashpot  C3.  Attached  to  the  main  mass  on  each 
side  are  two  coupled  beams  w  ith  tip  masses  m,  and  »k, 
stiffnesses  A',  and  A2,and  lengths  /,  and  /2.  In  the  analysis 


Fig.  I.  Schematic  diagram  of  an  aeroelastic  Structure 
and  Coordinate  System 


of  the  shown  system  only  the  symmetric  motions  of  the 
two  sides  of  the  model  are  considered.  Under  random 
excitation  the  system  response  will  be  described  by  the 
generalized  coordinates  <y , ,  q2 ,  and  q3  as  shown  in  the 
figure.  The  equations  of  motion  are  defived  by  applying 
Lagrange’s  equation 
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where  L=T—V. 

The  kinetic  energy  T  is  given  by  the  expression11 
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where  a  dot  denotes  differentiation  with  respect  to  time  t. 
Neglecting  the  gravitational  effects,  the  potential  energy 
V  is  eiven  bv 


l'=  1/2  (kp/i  +  k2qi+ktfi) 


(3) 


Substituting  for  T  and  V  in  equation  (1 ).  and  considering 
FU)  as  the  only  nonconservative  force  (damping  forces 
will  be  introduced  later)  results  in  the  equations  of  motion 
in  terms  of  the  ncndimensional  coordinates  <j,. 
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tf}  is  taken  as  the  root-mcan-squarc  of  the  mam  mass 
w  hen  all  other  parts  are  locked  under  forced  excitation. 
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o>3  is  taken  as  the  third  eigenvalue  of  the  system,  and 
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where  a  prime  denotes  differentiation  will)  respect  to  the 
dimensionless  time  r. 


III.  EIGENVALUES  OF  THE  SYSTEM 

The  system  eigenvalues  are  determined  from  the 
conservative  linear  part  of  the  equations  of  motion 

[>]{<?}  +  [>']{<?}  =  {°}  (6) 

The  characteristic  equation  of  (6)  is 

Dct|[/c]  — co3[/n]|  =  0  (7) 

where  a>  is  the  eigenvalue  of  the  mode  in  question. 
Expanding  the  determinant  gives  the  cubic  equation 


where  the  frequency  parameters  c o„=K,/m„.  are  the 
natural  frequencies  of  the  individual  components  of  the 
structure.  The  IMSL  (International  Mathematical  and 
Statistical  Library)  Subroutine  ZPOLR  (Zeros  of  a 
Polynomial  with  Real  Coefficients)  is  used  to  find  the 
roots  of  equation  (S)  numerically.  Fig.  2  shows  a  sample 
of  the  dependence  of  the  natural  frequency  ration/-  =  ui3/ 
(oi,  +m2)  on  the  ratios  M|  1/0/33  and  (o22/i:)33  for  beams 
length  ratio  l2  /7,  =  0.25,  and  mass  ratios  m2  ;'m,  =  0.5,  and 
imj/iii,  =  5.0.  Other  sets  of  curves  for  different  system 
parameters  are  obtained  and  reported  in  Ref  17  The 
importance  of  these  curves  is  to  define  the  critical  points 
where  the  structure  posscscs  internal  combination 
resonance  r=  1.0.  It  is  seen  that  the  most  critical  region  is 
located  for  the  curves  of  m22/w33  —  I  and  2.  For  the 
analysis  hereafter  the  following  parameters  will  be  used: 
/j//j  =  0.25,  w22/o)33  =  1.4. 

IV.  TRANSFORMATION  INTO  NORMAL 
COORDINATES 

Equations  (4)  include  linear  and  nonlinear  dynamic 
couplings.  The  linear  coupling  is  eliminated  In 
transforming  equations  (4)  into  normalized  cooidmates 
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by  using  the  transformation 
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where  [R]  is  the  modal  matrix  consisting  of  the 
normalized  eigenvectors, 
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the  elements  of  matrix  (10)  are  determined  by  using  the 
decomposition  method18  and  are  listed  in  Ref.  17. 

Rewriting  equations  (4)  in  the  matrix  form  and  using 
transformation  (9)  gives 

[>tt][R][Y")+[k][R]{Y)  =  :F)-{^  (11) 

Premultiplying  equation  (Ill  b\  the  transpose  of  the 
modal  matrix  results  in  diagonalizing  the  mass  and 
stiffness  matrices.  The  resulting  equations  unolvc 
nonlinear  coupling  and  have  the  form 
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V.  DYNAMIC  MOMENT  EQUATIONS 

The  response  coordinates  can  be  approximated  as  a 
Markov  vector  if  the  random  excitation  is  approximated 
as  a  zero  mean  physical  white  noise  U'(t)  having  the 
autocorrelation  function 

R„.(At)  =  E[\V(t)W{i  +  At)]  =  2D<5(At)  (14) 

where  2D  os  the  spectral  density  intensity  and  <5(  )  is  the 
Dirac  delat  funcation.  This  modelling  is  justified  as  long 
as  the  relevant  Wong-Zakai16  correction  term  is 
introduced.  The  non-lmear  functions  i //,  contain 
acceleration  terms  coupled  with  displacement 
coordinates  such  as  Y"Yj.  These  terms  arc  removed  from 
equations  (12)  by  successive  elimination  by  using 
MACSYMA  software.  Equations  (12)  take  the  new  form 

>7  +  2C,r,3  Y\  +  r/j  T,  =  /,U/(t)  +  Eg  AY,  Y)  (15) 

where  linear  viscous  damping  terms  have  ncen  introduced 
to  account  for  energy  dissip?  ion,  and 

cj2  =  (kjj/Af „■)(?;,  /in,),  /' , 3  =  w,/w3, 

/  =  hjM  e  =  qc}/l, 
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Introducing  the  transformation  into  the  Markov  state 
vector  X 

y].  y2<  r2,  y3,  >’1}  =  ;.y1..y: . a-6;  (i6> 

equations  (15)  may  be  written  in  the  standard  form  of 
Stratonovich  differential  equations 

d.Y,  =  f,(X.  r)dr  +  V  G,,(X.  r)dB,(r)  (17) 
/*  i 

where  the  white  noise  H'tr)  has  been  replaced  by  the 
formal  derivative  of  the  Brownian  motion  process  B( r). 
i  e  . 

(V(r)  =  <7  dfle  l  dr,  a2  ~2D 

Alternatively,  equations  (17)  may  m  turn  be 
transformed  into  the  Ito  type  equation 
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where  the  double  summation  expression  is  called  the 
Wong-Zakai  correction  term16. 

The  system  stochastic  Ito  equations  are 

dX",  =  X’2  dr,  dX’3  =  X'4  dr,  dXf5  =  A'6dT 
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The  evolution  of  the  response  probability  density 
function  is  described  by  the  Fokker-Planck  equation 
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where  p(X,r)  is  the  response  joint  probability  density 
function,  and  a,(X.  t )  and  6,;(X,  r)  arc  t he  first  and  second 


incremental  moments  evaluated  as  follows 
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The  coefficients  a,  and  b,j  are  evaluated  for  the  present 
system  with  the  aid  of  MACSYMA  program.  It  is  not 
possible  to  solve  the  resulting  Fokker  planck  equation 
even  for  the  stationary  case.  Instead,  one  may  generate  a 
general  first  order  differential  equation  describing  the 
evolution  of  response  moments  of  any  order.  This 
equation  is  obtained  by  multiplying  both  sides  of  the 
system  Fokker  Planck  equation  by  the  scalar  function 
<*>(X) 

<1>(X)  =  XV  X*--XVX*‘A'V  X  V  (22) 

and  integrating  by  parts  over  the  entire  state  space  —  co  < 
X<oo.  The  following  boundary  conditions  are  used 

p(X— »  — oo)  =  p(X— >oo)  =  0  (23) 

Due  to  space  limitation  the  system  moment  equation 
will  not  be  listed  in  this  paper.  The  reader  may  refer  to 
Ref.  1 7  for  more  details.  However,  the  general  form  of  the 
resulting  differential  equation  is 

>K'  =  F„(ni,,m2,. .  .,ms,ntN*i)  (24> 

where  N  =  £ f=  ,  kt. 

In  deriving  the  system  moment  differential  equation  the 
following  notation  is  adopted 


ac 

M:.  • 

J  -  c 


p(X, t)<3>(X)  dA' ,  dX2 


It  is  found  that  the  differential  equation  of  order  N 
contains  moment  terms  of  order  A'  and  A'  + 1 .  The  source 
of  this  infinite  hierarchy  is  the  system  nonlinear  functions 
t/q  in  equations  (12).  If  these  nonlinear  functions  are 
dropped  the  system  becomes  linear  and  the  response 
moment  equations  are  consistent.  In  the  present  study  the 
following  three  cases  will  be  examined: 

(i)  Linear  response  of  constant  coefficients  structure. 

(ii)  Linear  response  of  the  structure  with  random 
stillness. 

(iii)  Response  of  the  structure  with  autoparametric 
interaction  involving  the  internal  combination 
internal  resonance  023  =  02,  +02,. 


VI.  STRUCTURE  W  ITH  CONSTANT 
PARAMETERS 

The  equations  of  motion  for  this  case  are  obtained  from 
equations  (11)  by  excluding  the  nonlinear  functions  1 A,. 
The  resulting  equations  of  motion  arc 
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For  this  linear  case  the  response  moment  differential 
equations  are  consistent.  The  mean  squares  of  the 
stationary  response  is  obtained  in  the  closed  form 

E[Ytt  =  Df]  /( 2C ,  rf ,),  £[  >",J]  =  Df\  /( 2£ ,  r ,  3), 

£[  Y\]  =  Dfl  /(2f  2r|3),  £[  T^]  =  D/^  /(2£2r23), 

£[^]  E[  Y?]  =  Dfl/(2Z3) 

Before  presenting  the  linear  response  graphically,  it 
would  be  useful  to  recall  that  the  generalized  coordinates 
qt  were  nondimensionalized  with  respect  to  the  root-mean 
square  of  the  main  mass  response  when  the  coupled 
system  was  locked  under  forced  excitation.  The  value  of 
q3  can  be  estimated  from  the  single  degree  of  freedom 
equation  of  motion 

(w,  +  »i3+m3)q°3  +  C}qc3+k3q3  =  F(r)  (28) 

which  has  the  stationary  response 

£[^2]=£[?32]=A/2:3  (29) 

and  therefore 

=  (30) 

The  excitation  parameter  level  D/2Z3  is  chosen  so  that 
q°3  is  chosen  so  that  q\  is  unity  and  as  a  result  any 
deviation  from  unity  gives  a  measure  of  the  dynamic 
interaction  (linear  or  nonlinear)  with  other  mides.  For  the 
analysis  hereafter  the  excitation  level  will  be  chosen  such 
that 

D/2[3  =  1  (31) 

In  this  case  the  mean  square  response  (27)  is  reduced  to 


the  simple  form 
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Fig.  4.  Mean  square  response  of  normal  modes  for 

c1c2=o.o;,c3=o.o; 


Fig.  5.  Mean  square  response  of  generalized  coordinates 
for  Ci  =  C3  =  0.005,  £3  —  0.01 


The  linear  response  for  both  normalized  and  generalized 
coordinates  is  determined  for  various  damping  ratios. 
Figs  3  and  4  show  the  mean  square  responses  as  a 
function  of  the  frequency  ratio  r  for  two  sets  of  damping 
ratios.  It  is  seen  that  both  the  first  and  second  normal 
mode  mean  square  responses  decrease  faster  than  the 


Fig.  3.  Mean  square  response  of  normal  modes  for 

-o.o()5,  :3  =  o.oi 


Fig.  6.  Mean  square  response  of  generalized  coordinates 
for  c,  =C2  =0.01 ,  q3  —  0.0l 


third  mode  as  the  frequency  ratio  increases.  In  terms  of 
generalized  coordinates.  Figs  5  and  6  shows  that  t he  mean 
square  displacement  increases  while  the  two  beam 
displacements  decrease  with  the  frequency  ratio.  The  two 
sets  of  figures  show  the  well  known  control  damping  effect 
on  the  mean  square  responses. 


VII.  STRUCTURE  WITH  RANDOM  STIFFNESS 

The  equations  of  motion  of  this  case  are  obtained  by 
including  a  random  component  to  each  stifTness  in  the 
original  linear  equations  of  motion. 

The  equations  of  motion  take  the 


A' |  +  S,(r)  0  0 

0  Aj  +S2(/)  0 

0  0  A 


Aj  +  Sj(f) 


Introducing  the  same  dimensionless  parameters  listed 
in  Sections  II  and  III,  the  equations  of  motion  in  terms  of 
the  normal  coordinates  after  introducing  linear  damping 
are: 

.Vi+2;,)-13y;  +r?3y,  +s11W'1(t)+s12M/j(t) 

+  s13W3(x)]yl=flW(T) 

y'i  +  2 C2r23 V'i  +  rhy2  +  [s2 ,  W',  (t)  +  s22  W'2(t) 

+  s23IV3(T)]y2=/2lV(T) 

y'i  +  %  3y  3  +  )'3  +  [s3 1  IF,  (t  )  +  s32  w2  (t  ) 

+  5335F3(T)]y3=/3lV(t)  (34) 

where 

and  W,(t)  arc  zero  mean  white  noise  processes  with 
spectral  densities  20,.  Equations  (34)  constitute  a  set  of 
coupled  differential  equations.  The  response  mean 
squares  are  obtained  by  solving  the  stationary  moment 
equations.  The  analytical  solution  for  the  stationary 
response  is 

£[  )  l]  =  Df\  /{2,  |f  ?3  —  ~  f>2Sl2  ~  ^35l  3|  ’ 

E[Y?]  =  >-U£[r\] 

£[  Fj]  =  Dfl  /{2>2r23  —  D|521  —  D2s\2  —  D  }s\ 3] , 
E[Y22]  =  r\}E[Yl] 

£[  >  3]  =  D/j/{ —  D,s3 1  —  DjSj,  —  Dj553J, 

£[n:]  =  £[Fi]  (35) 

This  solution  indicates  that  the  system  may  be  unstable 
depending  on  the  values  of  D,.  The  fact  that  the  mean 
square  must  always  be  positive  provides  the  stability 
criteria  for  mean  squares  given  by  (35).  These  criteria  are 
obtained  by  keeping  the  denominators  of  (35)  always 
positive,  i.e., 

2£ir?3>  (0|S| ,  +D:s?2  +  D3s?3) 

2  Cir2)>(D,  M  i  +  D2s  2  2  +  D  j.s;  3 ) 

2> 3  >  i  +  D2s232  +  D (36) 

The  stability  boundaries  represented  by  conditions  (36) 
are  shown  in  Fig.  7  as  a  function  of  the  internal  resonance 
frequency  ratio  r.  For  simplicity  the  excitation  lev  els 
D,/2^,  of  the  random  stifTness  perturbations  arc  assumed 
to  be  equal  Samples  of  the  response  means  squares  as 
function  of  the  excitation  level  D: 22  are  shown  in  Fies  S 
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Fig.  7.  Mean  square  stability  boundary  of  the  structure 
with  random  stiffness,  for  ;,  =  0.01 


Fig.  S.  Mean  square  response  of  normal  modes  with 
random  stiffness  for  i  =  1.0,  f,  =  0.0I 


and  9  in  terms  of  normal  and  generalized  coordinates, 
respectively.  It  is  observed  that  the  response  of  tip  mass  of 
the  vertical  cantilever  is  the  main  source  of  instability. 

VIII.  AUTOPAR AMETR1C  INTERACTION 

In  this  case  the  influence  of  nonlinear  modal  coupling  on 
the  system  response  will  be  examined  by  including  the 
functions  in  the  analysis  These  functions  arc  only 
significant  if  the  structure  is  tuned  internally  such  that  the 
normal  mode  frequencies  have  a  linear  relationship.  For 
the  present  system  it  is  found  that  the  following  three 
internal  resonance  conditions  can  take  place' 5 : 

CJi  =  (U|  +  (1>2 

(d3=2  w,  and  <u3  =  2(u2  (37) 

The  random  response  of  the  system  will  be  cxmined 
under  the  first  internal  resonance  condition  As 
mentioned  in  Section  111  the  response  moment  equations 
involve  infinite  coupling  which  must  be  closed  in  order  to 
solve  for  the  response  statistics  It  is  known  that  the 
response  of  any  nonhneui  system  to  a  random  Gaussian 
excitation  will  he  non-Gaussian  The  deviation  of  the 
response  from  normality  depends  on  the  degree  of  the 
svstem  nonlincarnv  Genciallv.  closuie  schemes  me 
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Fig.  9.  Mean  square  response  of  generalized  coordinates 
for  same  conditions  of  FigsS 

classified  into  Gaussian  and  non-Gaussian16.  The 
Gaussian  schemes  are  useful  for  dynamic  systems  with 
weak  nonlinearity.  However,  in  certain  situations  the 
application  of  Gaussian  closures  may  lead  to  stochastic 
stability  boundaries  which  are  different  from  those 
derived  by  other  techniques  such  as  Stratonovich 
stochastic  averaging  or  non-Gaussian  closure 
approaches.  This  type  of  contradiction  has  been  reported 
for  nonlinear  systems  under  parametric  random 
excitations16.  For  two  degree -of-freedom  systems  the 
Gaussian  closure  scheme  yields  nonstationary  response 
while  non-Gaussian  closure  gives  strictly  stationary 
response  However,  the  main  response  characteristics  arc 
found  identical  as  predicted  by  both  methods. 

This  Section  examines  the  nonlinear  response  as 
obtained  by  using  a  Gaussian  closure  scheme  which  is 
based  on  the  properties  of  the  cumuiants.  For  the  present 
system  27  equations  for  the  first  and  second  order 
moments  will  be  generated.  The  moment  equations  are 
closed  by  setting  all  third  order  cumulants  to  zero,  i.e., 

;.3[.V,AVVJ  =  E[X,X,Xk]  - 1  E[X ,]£[XjA\] 

+  2£[.V,]£[.VJ]E[A't]  =  0  (38) 

where  the  number  over  summation  sign  refers  to  the 
number  of  terms  generated  in  the  form  of  the  indicated 
expression  without  allowing  permutation  of  indices. 
Relation  (38)  is  used  lo  obtain  expressions  for  the  third 
order  moments  in  terms  of  first  and  second  order 
moments. 

The  solution  of  the  closed  27  coupled  moment 
equations  is  obtained  numerically  by  using  the  1MSL 
DVERK  Subrouiine  (Runge-Kutia-Verncr  fifth  and  sixth 
numerical  integration  method)  Depending  on  the  value 
of  internal  detuning  parameter  r  the  system  response  may 
be  reduced  to  the  same  linear  response  of  section  VI  or 
may  become  quasi -stationary  which  deviates  significantly 
from  the  linear  solution  The  response  of  auloparamctric 
interaction  is  found  to  take  place  in  regions  of  internal 
resonance  ratio  slightly  deviated  from  the  exact  tuning 
r  =  I  The  dev  uition  may  be  attributed  to  the  contribution 
of  nonlmeaiiliCs  incurred  during  the  Gaussian  closure 
procedure  Surprisingly,  exact  internal  resonance  yields 


linear  response  characteristics  which  are  displayed  in  Fig. 
10.  It  is  seen  that  the  response  fluctuates  between  two 
limits  during  the  transient  period,  then  converges  to  a 
stationary  values  which  corresponds  exactly  to  the  linear 
solution  of  section  VI.  The  effect  of  different  initial 
conditions  is  examined  and  it  is  found  that  regardless  of 
the  initial  conditions  the  solution  reaches  the  same  steady 
state  value.  For  internal  resonance  ratio  r=  1.175,  Fig.  1 1 
shows  another  set  of  time  history  responses.  In  this  case 
the  response  mean  squares  do  not  achieve  a  stationary 
state.  During  the  transient  period  the  frequency  of  the 
third  mode  is  approximately  1 .17  times  the  sum  of  the  first 
two  mode  frequencies.  The  quasi-stationary  behaviour, 
although  present  for  all  three  modes,  is  most  prominent 
for  the  second  mode. 

To  further  illustrate  the  departure  of  the  nonlinear 
response  from  the  linear  one.  Figs  12-15  display  the 
dependence  of  the  normalized  mean  squares  on  the 
internal  resonance  for  various  system  parameters.  The 
mean  squares  arc  normalized  by  the  corresponding  linear 
solution.  The  subscript  G/L  refers  to  the  ratio  of  the 
nonlinear  Gaussian  solution  to  the  linear  response.  In  the 
regions  near  critical  internal  resonance  the  upper  and 
lower  envelopes  of  the  quasi-stationary  response  are 
plotted.  A  general  trend  is  observed  to  exist  in  all  figures. 
There  is  a  sharo  increase  in  the  displacement  mean  square 


Fig.  10.  Time  history  response  of  normal  coordinates  for 
£,  =  0.01 ,  c  =  0.025.  x-  013/(10,  +(Uj)=  I  0 
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Fig  II  Time  history  response  of  normal  com  Jinnies 
shonmg  autoparwnetne  interaction,  for  ',  =  (IHl .  1  - 
0.025.  r  —  (o3  ;(oj,  -t-  a>2 )  =  1 .175 
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Fig.  12.  Mean  square  resonse  of  normalized  coordinates 
as  function  of  internal  resonance  ratio  r,  for  £,  =  C2  = 
0.005.  i-L  —  0.01 ,  e  =  0.025 
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of  the  second  mode  associated  with  a  corresponding 
decrease  in  the  mean  square  of  the  third  normal  mode  and 
very  slight  drop  in  the  first  mode.  This  feature  is  similar  to 
a  great  extent  to  the  deterministic  nonlinear  absorbing 
effect  reported  by  Ibrahim  and  Woodall15.  Figs  12 and  13 
(  show  the  effect  of  damping  ratios  of  the  system  response. 
It  is  seen  that  nay  increase  in  damping  results  in 
narrowing  thr  region  of  autoparametric  interaction.  The 
nonlinear  coupling  parameter  c  has  a  direct  influence  on 
the  degree  of  the  response  deviation  from  the  linear 
solution  as  shown  in  Figs  12-15.  As  e  increases  from  0.025 
to  0.05  the  region  of  autoparametric  interaction  becomes 
i  more  wider. 

IX.  CONCLUSIONS 

The  linear  and  nonlinear  modal  interactions  of  a  three- 
degree-of-freedom  structure  subjected  to  random 
excitation  is  examined.  For  the  linear  modelling  the 
response  is  determined  for  two  cases  of  structure 
parameters.  The  first  case  is  when  the  parameters  are 
constant  coefficients.  The  mean  square  response  of  this 
case  is  obtained  in  terms  of  the  excitation  spectral  density 
and  the  internal  detuning  parameter.  The  second  case 
involves  random  parametric  excitations  in  the  stiffness 
matrix.  These  excitations  result  in  modal  parametric 
coupling  of  the  normal  coordinates.  The  mean  square 
responses  are  governed  by  the  spectral  densities  of 
parametric  excitations  which  also  result  in  the  conditions 
of  mean  square  stability.  The  results  of  the  first  case  are 
used  as  a  reference  to  measure  the  effects  of  nonlinear 
inertia  coupling  of  normal  modes  on  the  mean  square 
response  of  the  system  in  the  neighbourhood  of 
combination  internal  resonance.  It  is  found  that  the 
critical  internal  resonance  occurs  at  a  value  close  to  r  = 
1.175  which  is  deviated  from  the  exact  value  r=  1.  The 
nonlinear  modal  interaction  results  in  an  increase  of  the 
second  normal  mode  mean  square  response  and  in  an 
associated  decrease  of  the  first  and  third  normal  modes. 

ACKNOWLEDGEMENT 

This  research  is  supported  by  a  grant  from  the  Air  Force 
Office  of  Scientific  Research  under  grant  no.  AFOSR-S5- 
0008  Dr  A.  Amos  is  the  Proaram  director. 


REFERENCES 

1  Navfeh.  A  H  and  Mook.  D  T.  Nonlinear  Oscillations,  John 
Wiley  and  Sons,  New  York,  1979 

7  Schmidt.  G  and  Tondl,  A  Nonlinear  Vibration ,  Cambridge 
University  Press,  Cambridge.  1985 

3  Breilbach.  E.  J.  Effect  of  structural  nonltnearitics  on  aircraft 
vibration  and  flutter.  Advisroy  Group  for  Aerospace  Research 
and  Development  of  NATO,  AGARD  Report  No  665. 
September  1977 

4  Breilbach.  E.  J  Flutter  analysis  of  an  airplane  with  multiple 
structural  nonlinearities  in  the  control  system,  NASA  Technical 
Paper  1620.  1980 

5  De  Ferrari.  G..  Chcsia,  L..  Sensburg.  O.  and  Lotae.  A.  Effects  of 
nonlinearity  on  wing-store  flutter.  Advisory  Group  for 
Aerospace  Research  and  Development  of  NATO,  AGARD-R- 
687.  I9S0 

6  Pel.  'ubcl.  R.  P  .  Jr.  Haller.  R  L  and  Bolding.  R.  M  F- 16  Flutter 
suppression  system  investigation.  Proceedings  of  the  /II A  4, 
ASM  £fASCE:  AHS  7 /si  Structures.  Structural  Dynamics  ami 
Materials  Conference.  Seattle,  Washington.  1980.  620-634 

7  Reed,  W.  H  111,  Foughner,  J.  T..  Jr  and  Runyan.  H  L. 
Decoupler  pvlon:  A  simple  effective  wing/sotre  suppressor,  J. 
Aircraft.  19S0,  17(3).  206-211 

8  Desmarais.  R.  N.  and  Reed,  W.  H  .  Ill  Wing/store  flutter  with 
non-linear  pylon  stoffness,  J .  Aircraft.  1981,  18(1  J ),  984-987 

9  Bolotin.  V.V.  The  Dynamic  Stability  of  Elastic  Sy stems,  Holden- 
Day.  Inc..  SanFrancisco.  1964 

10  Barr,  A.  D.  S  Some  developments  in  parametric  stability  and 
nonlinear  vibration.  Proceedings  of  International  Conference  of 
Recent  Advances  in  Structural  Dynamics.  Southampton, 
England.  1980,  545-568 

1 1  Minorskv,  N.  Nonlinear  Oscillations.  Van  Nostrand,  New  York, 
1962 

12  Barr,  A.  D.  S  and  Ashworth,  R  P.  Parametric  and  nonlinear 
mode  interaction  behavior  in  the  dynamics  of  structures, 
AFOSR  Final  Report,  Department  of  Mechanical  Engineering. 
University  of  Dundee,  December  1977 

13  Haddow.A  G  .  Barr,  A.  D.  S.  and  , Mook,  D  T.  Theoretical  and 
experimental  study  of  modal  interaction  in  a  two-degree-of- 
freedom  structure.  J.  Sound  and  Vibration,  1984,  97(3),  451-473 

14  Ibrahim,  R.  A  .  Woodall.  T  D  and  Heo,  H  Modal  analysis  of 
structural  systems  involving  nonlinear  coupling.  Shock  and 
Vibration  Bulletin.  1984.  54(1),  19-27 

15  Ibrahim.  R  A.  and  Woodall,  T.  D  Linear  and  nonlinear  modal 
analysts  of  aeroclastic  structural  systems.  Computers  and 
Structures,  19S6,  22(4),  699-707 

16  Ibrahim.  R.  A  Parametric  Random  Vibrations,  John  Wiles  and 
Sons,  New-  York.  1985 

17  Hedayati,  Z  Random  modtl  interaction  of  a  nonlinear 
aeroelastic  structure,  MS  Thesis,  Texas  Technical  University. 
August  1985 

IS  Craig.  R  R  Structural  Dynamics  -  An  Introduction  lo  Computer 
Methods.  John  Wiley  and  Sons.  New  York.  1981 


Stochastic 

modal  interaction: 

R. 

A. 

Ibrahim 

and 

Z 

Hedayati 

Stochastic 

modal  interaction 

R. 

A. 

Ibrahim 

and 

7 

Hedayati 

Stochastic 

modal  interact  ion 

R 

A. 

Ibrahim 

and 

Z 

Hedayati 

Slot  hast  tc 

modal  interact  ton 

R 

A. 

Ibrahim 

and 

2. 

Hedayati 

Stochastic 

modal  interaction 

R 

A 

Ibrahim 

and 

2. 

Hedayati 

Stochastic 

modal  interact  ion 

R 

A. 

Ihralum 

and 

2 

Hedayati 

Stochastic 

modal  interact  ton 

R 

A. 

Ibrahim 

and 

2 

Hedayati 

Stochastic 

modal  interaction 

R 

A. 

1  hr  alum 

and 

2 

Hedavati 

Stochastu  modal  interim  ion  R  4  Ibrahim  anil  2  Hedayati 


EXPERIMENTAL  INVESTIGATION  OF  STRUCTURAL  AUTOPARAMETRIC 


INTERACTION  UNDER  RANDOM  EXCITATION 


R.  A.  Ibrahim*  D.  C.  Sullivan** 
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Lubbock,  Texas  79409 


ABSTRACT 

The  paper  presents  the  results  of  an  experimental 
investigation  of  random  excitation  of  a  nonlinear 
tvo-degree-of -freedom  structural  model.  The  model 
normal  mode  frequencies  are  adjusted  to  have  the 
ratio  of  2  to  1.  This  ratio  meets  the  condition  of 
internal  resonance  of  the  analytical  model.  When 
the  first  normal  mode  is  externally  excited  by  a 
band  limited  random  excitation,  the  system  mean 
square  response  is  found  to  be  linearly  propor¬ 
tional  to  the  excitation  spectral  density  up  to  a 
certain  level  above  which  the  two  normal  modes 
exhibit  discontinuity  governed  mainly  by  the  inter¬ 
nal  detuning  parameter  and  the  system  damping 
ratio.  The  results  are  completely  different  when 
the  second  normal  mode  is  externally  excited.  For 
small  levels  of  excitation  spectral  density  the 
response  is  dominated  by  the  second  normal  mode. 
For  higher  levels  of  excitation  spectral  density 
the  first  normal  mode  attends  and  Interacts  with 
the  second  normal  mode  in  a  form  of  energy 
exchange.  A  number  of  deviations  from  theoretical 
results  are  observed  and  discussed. 


I.  IFTBPPVCIIOP 

The  last  two  decades  have  witnessed  an  Increasing 
interest  in  the  study  of  dynamic  behavior  of  non¬ 
linear  systems  under  deterministic  and  random  exci¬ 
tations.  Under  certain  conditions  these  systems  may 
experience  complex  response  characteristics  such  as 
jump  phenomenon,  limit  cycles,  internal  resonance, 
saturation  phenomenon,  and  chaotic  motion.  These 
nonlinear  phenomena  have  been  predicted  theoreti¬ 
cally  ~Jand  observed  experimentally4*6  under  harmo¬ 
nic  excitations.  However,  most  of  the  predicted 
random  response  characteristics,  including  response 
stochastic  stability  and  statistics,^’®  have  not 
been  verified  experimentally.  Very  few  experimental 
investigations  of  random  vibration  of  nonlinear 
systems  have  been  reported  in  the  literature.^  The 
lack  of  experimental  verifications  may  be  due  to 
several  reasons.  These  include  difficulties  in 
generating  the  same  properties  of  the  random  exci¬ 
tation  as  represented  theoretically,  and  the  limi¬ 
tations  of  experimental  equipment.  Recently, 
Bolotin"  discussed  a  number  of  experimental  diffi¬ 
culties  encountered  in  experimental  measurements  of 
stochastic  stability  of  parametric  excited  systems. 

In  deterministic  nonlinear  vibrations;  the  ampli¬ 
tude  Jump,  limit  cycles,  and  parametric  instability 
are  common  features  of  nonlinear  single-  and  multl- 
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degree -of- freedom  systems.  Parametric  instability 
takes  place  when  the  external  excitation  appears  as 
a  coefficient  in  the  homogeneous  part  of  the  equa¬ 
tion  of  motion.  It  occurs  when  the  excitation 
frequency  is  twice  (or  multiple)  of  the  system 
natural  frequency.  Internal  resonance  and  satura¬ 
tion  phenomenon  may  occur  only  in  nonlinear  systems 
with  more  than  one  degree-of-freedom.  Internal 
resonance  implies  the  existence  of  a  linear  rela¬ 
tionship  between  the  system  natural  frequencies  and 
causes  nonlinear  normal  mode  interaction  in  the 
form  of  energy  exchange.  Under  external  harmonic 
excitation,  the  mode  which  is  directly  excited, 
exhibits  in  the  beginning,  the  same  features  of  a 
single-degree-of  freedom  system  response  and  all 
other  modes  remain  dormant.  As  the  excitation 

amplitude  reaches  a  certain  critical  level,  the 
other  modes  become  unstable  and  the  originally 

excited  mode  reaches  an  upper  bound.  In  this  case, 
the  mode  is  said  to  be  saturated  and  energy  is 
transferred  into  other  modes.  This  interesting 
phenomenon  takes  place  only  in  systems  with  quadra¬ 
tic  nonlinear  coupling  which  results  in  a  third 
order  internal  resonance. 

Under  deterministic  unsteady  aerodynamic  forces, 
most  nonlinear  characteristics  can  be  predicted  by 
one  of  the  standard  techniques  of  nonlinear  diffe¬ 
rential  equations.  However,  aerospace  structures 
are  usually  subjected  to  turbulent  air  flow,  and 
the  aeroelastician  is  confronted  with  aerodynamic 
loads  which  are  random  in  nature.  These  loads  vary 
in  a  highly  irregular  fashion  and  can  be  described 
in  terms  of  statistical  quantities  such  as  means, 
mean  squares,  autocorrelation  functions  and  spec¬ 
tral  density  functions.  Ibrahim  and  Roberts  ’ 
and  Ibrahim  and  Heo l2’  ^considered  nonlinear  two 
degree-of-freedom  structural  systems  and  applied 
Gaussian  and  non-Gaussian  closure  techniques  to 
predict  the  response  statistics  and  response  sto¬ 
chastic  stability.  These  studies  revealed  that  a 
system  with  internal  resonance  may  experience  non¬ 
linear  characteristics  such  as  autoparametrlc  in¬ 
teraction.  Roberts  14  conducted  a  series  of  expe¬ 
rimental  tests  to  measure  the  mean  square  stability 
boundaries  of  a  unlmodal  response  of  a  coupled  two - 
degree-of-freedom  system.  Roberts  reported  a  num¬ 
ber  of  difficulties  in  measuring  the  stability 
boundaries.  Based  on  the  authors  experience  and 
other  Investigators  work,  it  is  understood  that 
experimental  investigation  of  non'' tnear  random 
vibration  is  not  a  simple  task  and  requires  careful 
planning  and  advanced  equipment  preparations. 

The  purpose  of  the  present  paper  is  to  report  the 
results  of  an  experimental  investigation  to  mea¬ 
sure  the  response  mean  squares  of  a  nonlinear  two 
degree-of-freedom  structural  model  under  band 
limited  random  excitation.  The  same  model  was  ana¬ 
lytically  examined  by  Haddov,  at  al.  under  harmo- 
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nlc  excitation,  and  by  Ibrahim  and  Heo1^*^  under 
wide  band  random  excitation.  Agreements  and  disag¬ 
reements  with  theoretical  predictions  will  be  dis¬ 
cussed  together  with  recomoendations  for  future 
experimental  work. 


The  random  response  of  a  two  degree-of-freedom 
elastic  structure  has  been  determined  analytically 
In  references  [12,13].  The  analytical  model  shown 
In  fig.  (1)  consists  of  two  beams  with  end  masses. 
Under  vertical  support  motion  £(t)  the  response  of 
the  two  beams  Is  mainly  governed  by  linear  dynamic 
and  parametric  couplings.  However,  if  the  system 
is  designed  such  that  the  first  two  normal  mode 
frequencies  u  and  u  satisfy  the  internal  reso¬ 
nance  condition  id-  *  2u>^  ,  the  nonlinear  inertia 
forces  become  dominant  and  the  system  dynamic  res¬ 
ponse  deviates  from  the  linear  response.  In  terms 
of  the  non-dimensional  normal  coordinates  Y  the 
system  equations  of  motion  are: 


[I ] (Y" )  +  [ C] (Y*  >  +  ( r2 ] ( Y ) 


£  "(t) (a)  +  e£-(T)[b](Y)  +  £<  *  )  (1) 

where  a  prime  denotes  differentiation  with  respect 
to  the  nondimens ional  time  parameter  T  -  wt .  and  the 
coordinates  Y  are  related  to  the  dimensional  normal 
coordinates  y  by  the  relation  (Y^  .Y^  )  “  • 

q°  is  taken  as  the  response  root  mean  square  of  the 
system  when  the  length  of  the  vertical  beam  shrinks 
to  zero,  l.e.  the  response  root  mean  square  of  the 
main  beam  with  end  mass  (m  +m  ) .  The  elements  of 
the  vector  (a)  and  matrix  [bl  are  constants  depen¬ 
ding  on  the  system  properties.  The  small  parameter 
e-  q° /l  .  The  matrix  [ir]  is  diagonal  with  elements 
1  ani  (A  /w  )2 •  The  vector  (  y  )  contains  all  quad¬ 
ratic  noftliitear  terms  which  encompasses  two  groups: 
nonlinear  terms  of  the  same  mode  and  autoparametric 
terms  of  the  type  Y  Y"  It  is  the  autoparametric 
coupling  which  gives  tlse  to  the  internal  resonance 
condition  r  -  ”  2' 

The  random  acceleration  Z\t)  was  assumed  to  be 
Gaussian  wide  band  process  with  zero  mean  and  a 
smooth  spectral  density  2D  up  to  some  frequency 
higher  than  any  characteristic  frequency  of  the 
system.  The  acceleration  terms  in  the  nonlinear 
functions  iK  were  removed  by  successive  elimina¬ 
tion  and  the  system  equations  of  motion  was  tran¬ 
sformed  into  a  Markov  vector  via  the  coordinates 
transformation 
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A  set  of  first  order  differential  equations  of  the 
response  statistical  moments  were  generated  by 
using  the  Fokker-Planck  equation  approach.'  These 
equations  were  found  to  be  coupled  through  higher 
order  moments  and  were  closed  via  two  approaches: 
Gaussian  and  non-Gaussian  closures.  These  closure 
techniques  are  based  on  the  cumulant  properties. 
The  Gaussian  closure  is  established  by  equating  all 
cuaulants  A  of  order  greater  than  two  to  zero,  l.e. 


ql 


Fig.(l)  Model  of  coupled  two  beams 


Fig.  (2)  Gaussian  closure  solution  tor 
various  values  of  nonlinear 
coupling  parameter  e 


This  approach  resulted  in  fourteen  coupled  differe¬ 
ntial  equations  for  first  and  second  order  moments 
of  the  response  coordinates.  The  numerical  integ¬ 
ration  of  these  equations  revealed  that  the  respon- 
te  mean  squares  fluctuate  between  two  limits.  This 
fluctuation  means  that  the  response  does  not 
achieve  a  stationary  state.  The  autoparametric 
interaction  took  place  In  the  neighborhood  of 
internal  resonance  and  was  manifested  by  an  energy 
exchange  between  the  mean  squares  of  the  two  normal 
modes.  Figure  (2),  taken  from  reference  12,  shows 
a  sample  of  the  mean  square  response  of  the  system 
normal  modes  against  the  internal  detuning 
parameter. 

The  second  method  takes  into  account  the  effect  of 
the  response  non-normality.  As  a  first  order  non- 
Gaussian  approximation  all  cumulants  of  order 
greater  than  four  were  equated  to  zero,  l.e. 
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>ig.  (3)  Non-Gaussian  closure  solution  ior 

various  values  oi  nonlinear  coupling 
parameter  c 


IXX-  x!°>-0'  >-5.^  «> 

This  approach  resulted  in  69  first  order  differen¬ 
tial  equations,  in  the  first  through  the  fourth 
order  moments,  which  were  solved  numerically.  The 
solution  reaches  a  stationary  state  after  a  tran¬ 
sient  period  and  exhibits  the  same  nonlinear 
interaction  as  predicted  by  the  Gaussian  closure 
solution.  Figure  (3)  shows  the  stationary  mean 
square  response  of  the  normal  coordinates  against 
the  Internal  detuning  parameter. 

Although  the  two  approaches  yielded  common  features 
to  those  predicted  by  deterministic  theory^  of 
nonlinear  vibration  such  as  autoparametrlc  sup¬ 
pression  effect,  the  random  analysis  did  not  verify 


the  existence  of  saturation  phenomenon.  The  satu¬ 
ration  phenomenon  ^  is  a  well  known  feature  for 
multi -degree -of- freedom  systems  involving  quadratic 
nonlinear  coupling  subjected  to  harmonic  forced 
excitation. 

It  is  well  known  that  the  predicted  results  are 
approximate  and  their  validity  has  not  been 
examined.  The  next  section  reports  the  measured 
results  of  a  series  of  experimental  tests  of  the 
same  model  under  band  llmltad  random  excitation. 


III.  EXPERIMENTAL  INVESTIGATION 


III ■ 1  Experimental  Model  and  Equipment 

The  model  is  similar  to  a  great  extent  to  the 
experimental  model  used  by  Haddow,  et  al .  It  con¬ 
sists  of  a  horizontal  beam  of  cross  section  of 
0 . lll"xl . 0" ,  length  7.5",  and  carries  a  tip  mass  of 
0.015  slug.  The  tip  mass  has  a  provision  for 
clamping  the  vertical  beam  which  has  cross  section 
0.054"xl.0".  The  length  of  the  vertical  beam  can 
be  adjusted  by  changing  the  location  of  its  top 
mass  (0.0127  slug).  The  deflections  of  the  two 
beams  are  measured  by  strain  gages  fixed  at  the 
root  of  each  beam.  Two  gages  are  mounted  on  the 
horizontal  beam  in  a  two  arm  bridge .  Four  gages 
are  mounted  on  the  vertical  beam  in  a  four  arm 
bridge.  The  fixed  end  of  the  horizontal  beam  is 
clamped  by  a  fixture  which  is  bolted  on  the  top  of 
the  shaker  armature.  The  shaker  is  a  Caldyne  model 
A88  of  thrust  100  lb  and  provides  1"  peak-to-peak 
stroke.  The  shaker  is  powered  by  a  Ling  Electro¬ 
nics  Model  RA-250  power  supply  and  receives  a  ran¬ 
dom  signal  through  a  GenRad  Type  1381  Random  Noise 
generator.  The  random  signal  is  filtered  to  a 
desired  band  width  with  a  Krohn-Hite  Model  3)63 
Variable  Electric  Filter.  The  filtered  signal  is 
amplified  via  a  Calex  Model  176  Instrument 
Amplifier.  Figure  (6)  shows  a  schematic  diagram  of 
the  instrumentation  used  in  this  investigation. 
The  acceleration  of  the  shaker  platform  is  measured 
by  a  FCB  Plezotronlc  Model  302A02  shock  accelerome¬ 
ter.  The  accelerometer  is  powered  by  a  PCB 
Piezotronic  Model  680C06  power  unit. 

The  first  two  normal  mode  frequencies  of  the  system 
are  determined  theoretically  and  measured  experime¬ 
ntally  as  a  function  of  the  beams  length  ratio  £  A. 
as  shown  in  fig.  (5).  This  figure  shows  that  xheJ 
internal  resonance  w  /wj  -  2  is  obtained  in  two 
locations  of  the  length  ratio.  At  these  mass 
locations  the  normal  mode  frequencies  are: 
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C/4  -  0.485: 

f  -  9.1  Hz. 

f  -  18.2  Hz 

(5a) 

2  1 

1 

2 

t2/i1  -  0.707: 

fj  -  7.45  Hz,  f2  -  14.9  Hz 

(5b) 

The  analog  signals  of  Che  exclCation  and  responses 
are  read  and  converted  into  binary  numbers  using  a 
Data  Translation  Model  DT-3752  Intelligent  Analog 
Peripheral  (IAP).  This  ZAP  is  capable  of  reading 
either  8  channels  (+10v)  or  16  channels  (0-10v)  of 
input.  It  can  also  read  and  convert  analog  signals 
at  up  to  40k  points  per  second.  This  unit  is 
mounted  in  an  expansion  slot  of  an  IBM  System  9001 
Benchtop  Computer.  The  control  and  programming  of 
the  Analog/Dlgltal  (A/D)  system  are  accomplished 
through  the  software  controlled  registers  and  field 
selectable  (hardware)  options.  The  software  con¬ 
trolled  registers  are  the  control  registers,  status 
register,  and  gain/channel  register.  The  control 
register  controls  the  operation  and  mode  of  the  A/D 
system.  The  modes  which  are  used  in  this  investi¬ 
gation  are  direct  memory  transfer  and  Increment 
mode  operation.  Direct  memory  transfer  places 
converted  data  directly  into  the  memory  of  the 
computer.  The  increment  mode  allows  the  A/D  to 
Increment  the  input  channel  number  automatically 
before  each  A/D  conversion.  This  allows  data  to  be 
taken  from  sequential  channels  without  requiring  a 
program  to  specify  each  channel.  The  status  regis¬ 
ter  reports  the  complete  status  of  the  A/D  system 
during  the  operation.  The  gain/channel  register 
selects  the  desired  channels  from  which  the  data  is 
to  be  taken  and  sets  a  programmable  gain  for  all 
input  signals.  This  gain  is  set  to  one  for  all 
tests.  The  computer  controls  the  DT-3752  through  a 
Fortran  program.  Analog  signals  are  converted  for 
a  specified  amount  of  time  or  until  the  computer 
memory  is  full.  When  the  computer  has  completed 
collecting  data,  the  data  is  transferred  to  a  flop¬ 
py  disk  for  future  processing. 

The  data  processing  is  performed  at  equally  spaced 
Intervals.  The  problem  of  determining  this  time 
interval  is  well  discussed  in  Bendat  and  Piersol'5 
Cenerally,  if  sampling  is  prepared  at  points  which 
are  too  close  together,  it  will  yield  correlated 
and  redundant  data.  This  will  unnecessarily  in¬ 
crease  the  labor  and  cost  of  calculations.  Sam¬ 
pling  at  points  which  are  too  far  will  lead  to  the 
problem  of  aliasing  The  aliasing  is  mainly  a 
confusion  between  the  low  and  high  frequency  compo¬ 
nents  in  the  original  data.  In  order  to  eliminate 
the  problem  of  aliasing,  a  sampling  rate  should  be 
chosen  to  be  at  least  two  time  the  maximum  frequen¬ 
cy  that  the  model  will  experience.  In  order  to  get 
a  good  sample  data,  a  sampling  rate  is  chosen  which 
is  roughly  eight  times  the  maximum  frequency.  In 
the  present  investigation,  the  sampling  rate  is 
chosen  to  be  80  Hz  per  channel  for  the  first  mode 
excitation  and  160  Hz  per  channel  for  the  second 
mode  end  wide  band  excitation.  Data  processing 
involves  another  problem  known  as  quantization 
which  is  the  conversion  of  data  values  at  the 
sampling  points  into  dlgitsl  form.  The  infinite 
number  of  values  of  the  continuous  analog  signal 
must  be  approximated  by  a  fixed  set  of  digital 
levels.  A  choice  between  two  consecutive  levels 
will  be  required  becsuse  the  scsle  is  finite.  The 
accuracy  of  the  approximating  process  is  a  function 
of  the  available  levels  which  is  dependent  upon  the 
analog  to  digital  converter  resolution.  The  accu¬ 
racy  of  the  DT-3752  is  the  value  of  the  least 
significant  bit  which  corresponds  to  a  voltage  of 


+0.0049v.  This  resolution  Is  analogous  to  a  def¬ 
lection  of  the  horizontal  beam  bean  of  +0. 00073- in 
and  the  vertical  of  +0.00097 -in  and  an  acceleration 
of  +0. 00044- g  for  the  excitation. 

The  experimental  model  is  tested  under  vsrious 
levels  of  excitation  spectral  density.  This  is 
achieved  by  keeping  the  input  signal  level  constant 
(Master  Gain  on  Ling  Amplifier)  for  the  range  of 
internal  detuning  of  the  model.  The  level  of 
amplification  is  adjusted  to  five  levels  for  tes¬ 
ting  of  both  the  first  and  second  normal  frequency 
bandwldths.  Another  series  of  tests  are  conducted 
for  excitation  spectral  density  that  covers  both 
normal  mode  frequencies. 

III. 2  Experimental  Results 

The  experimental  results  Include  sample  records  of 
time  history  responses  and  the  mean  square  respon¬ 
ses  in  terms  of  generalized  coordinates  and  normal 
coordinates.  The  mean  square  response  will  be  rep¬ 
resented  against  the  internal  detuning  parameter  r 
and  the  excitation  spectral  density  level.  The 
bandwidth  of  the  random  excitation  depends  on  the 
mode  under  investigation. 


III.  2.1  First  Mode  Excitation 

The  first  mode  is  excited  by  a  limited  bandwidth 
random  excitation  of  bandwidth  5Hz  and  a  central 
frequency  very  close  to  the  first  normal  mode  natu¬ 
ral  frequency.  The  frequency  content  of  this  ran¬ 
dom  process  is  selected  such  that  it  does  not 
excite  any  higher  structural  modes.  For  the  five 
levels  of  excitation  spectral  density,  the  system 
response  is  governed  mainly  by  the  first  mode  which 
does  not  show  any  nonlinear  coupling.  Figure  (6) 
shows  a  sample  of  the  time  history  response  under 
excitation  spectral  level  So-  0.0142  (g^/Hz)  when 
the  model  is  internally  tuned  to  the  resonance 
condition  /w.  “2.0.  It  is  seen  that  the  response 
is  characterized  by  a  narrow  band  random  process  of 
frequency  close  to  the  first  normal  mode  -  7.5  Hz. 


Fig.  (6)  Time  history  response  of  first  normal  mode 
excitation,  level  V,  Sq-  0.0142  g  ^/Hz 


Figure  (7a)  shows  the  mean  square  response  of  the 
generalized  coordinates  for  the  same  excitation 
spectral  density  level  of  fig.  (6).  The  empty 
points  are  measured  when  the  mass  of  the  vertical 
beam  moves  upward  while  the  full  points  are 
obtained  when  the  mass  moves  downward.  Both  groups 
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Fig  (7a)  Mean  square  responses  of  generalized 
coordinates  under  first  normal  mode 
excitation. 


Fig.  (7b)  Mean  square  response  of  normal 

coordinates  under  first  normal  mode 
excitation. 

Excitation  level  V:  S 0 -  0.0142  g  2/Hz .  empty  points 
correspond  nigber  position  of  the  upper  mass  while 
full  points  correspond  lower  position. 


are  measured  In  the  neighborhood  of  the  system 
Internal  resonance.  The  group  of  full  points  indi¬ 
cates  that  the  mean  square  of  the  horizontal  beam 
Increases  while  the  mean  square  of  the  vertical 
beam  decreases  as  the  normal  mode  frequency  ratio 
u  /w.  increases.  This  implies  that  the  model 
behaves  like  a  single  degree -of -freedom  system  for 
u2/ujj  »  2.  For  the  second  group  of  results  (empty 
points)  the  mean  square  response  of  the  vertical 
beam  Increases  and  the  mean  square  of  the  horizon¬ 
tal  beam  decreases.  This  feature  is  belonging  to 
the  characteristics  of  linear  vibration  absorbers 
due  to  inertia  coupling.  The  corresponding  respon¬ 
se  curves  In  normal  coordinates  are  shown  in  fig. 
(7b).  The  square  points  (empty  or  full)  are  belon¬ 
ging  to  the  first  normal  mode  which  obviously  pre¬ 
dominates  the  response.  It  is  also  seen  that  as 
the  vertical  mass  moves  downward,  the  model  starts 
to  behave  like  a  linear  single  degree-of-freedoo 
system  whose  mean  square  is  given  by  the  relation¬ 
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Fig  (8a)  Relationship  between  mean  square 

responses  of  generalized  coordinates 
and  the  excitation  spectral  density  for 
various  values  of  internal  detuning. 
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Fig  (8b)  Relationship  between  mean  square 

responses  of  normal  coordinates  and 
the  excitation  spectral  density  for 
various  values  of  internal  detuning. 


where  m,  u>  and  t  are  the  mass,  natural  frequency, 
and  damping  ratio  of  the  system,  respectively.  2D 


(Measurements  are  taken  for  lower  position  of  the 
upper  mass) . 


1*  the  excitation  spectral  density  of  a  wide  band 
random  excitation. 

It  is  clear  that  the  trend  of  the  full  square 
points  agrees  with  the  linear  solution  (6)  that  the 
mean  square  response  is  inversely  proportional  to 
the  cube  of  the  first  normal  mode  frequency. 

In  order  to  provide  more  insight  to  the  system 
response  statistics ,  the  mean  square  response  is 
plotted  against  the  excitation  spectral  density 
level  as  shown  in  fig.  (8a)  for  various  values  of 
internal  detuning.  It  is  seen  that  the  mean 
squares  of  the  two  beams  Increase  with  the  excita¬ 
tion  spectral  density  up  to  a  certain  level  above 
which  the  curves  are  discontinuous .  The  degree  of 
discontinuity  depends  on  the  internal  detuning. 
Any  deviation  from  the  exact  internal  detuning 
results  in  a  strong  discontinuity.  This  disconti¬ 
nuity  means  that  the  system  is  unstable  in  the  mean 
square  sense.  Similar  features  were  reported  in 
the  deterministic  response  of  the  same  system  by 
Haddow,  et  al.5  The  location  of  discontinuity  is 
strongly  dependent  on  the  values  of  damping  ratios 
and  the  internal  detuning  of  the  structure.  Figure 
(8b)  shows  the  mean  square  response  of  the  normal 
coordinates  against  the  excitation  spectral  densi¬ 
ty.  The  curves  have  the  same  trend  of  fie. 


The  second  normal  mode  is  excited  by  a  limited  band 
random  excitation  of  bandwidth  5  Hz  and  central 
frequency  very  close  to  the  second  normal  mode 
frequency.  Five  levels  of  excitation  spectral 
density  ranging  from  0.001  g2/Hz  to  0.022  g2/Hz  are 
selected.  A  general  feature  of  the  time  history 
response  records  is  that  both  amplitudes  q  and  q 
increase  with  the  levels  of  excitation  as1  in  thi 
first  mode  excitation.  The  records  also  show  that 
for  all  selected  beam  length  ratios  and  for  all 
levels  of  excitation  spectral  density,  the  vertical 
beam  amplitude  q  is  always  greater  than  the  hori¬ 
zontal  beam  amplitude  a  .  Another  observation  is 
that  when  the  excitation  level  is  held  constant  the 
amplitudes  q,  and  q2  Increase  slightly  as  the  beam 
length  ratio  Increases.  For  small  levels  of  exci¬ 
tation  spectral  density,  the  second  normal  mode  is 
observed  to  have  no  interaction  with  the  first 
mode.  However,  above  a  certain  level  of  excitation 
spectral  density,  it  is  found  that  the  first  mode 
appears  for  a  certain  period  of  time  and  then 
disappears  as  the  second  mode  takes  over,  and  so  on 
as  shown  in  fig.  (9a).  This  nonlinear  Interaction 
of  the  two  normal  modes  is  more  clarified 
in  fig.  (9b).  Under  harmonic  excitation,  Haddow, 
et  al.^  reported  similar  energy  exchange  between 
the  two  modes.  Furthermore,  it  was  shown  that  the 
directly  excited  mode  becomes  saturated  and  energy 
is  transferred  to  the  first  mode.  In  the  present 
investigation,  the  energy  transfer  takes  place  not 
only  under  high  levels  of  excitation  spectral  den¬ 
sity  but  also  when  the  the  internal  resonance  is 
approaching  the  value  2  as  vertical  beam  length  is 
increasing. 


The  mean  square  responses  of  the  generalized  and 
normal  coordinates  are  plotted  against  the  Internal 
detuning  parameter  r  in  figs.  (10a)  end  (10b), 
respectively.  The  suppression  effect  of  the 
excited  mode  takes  place  only  whan  the  vertical 
■•»»  is  moved  downward  as  shown  in  fig.  (10b)  by 
the  full  triangular  points.  The  second  node  mean 
square  (empty  triangular  points)  Increases  with  a 
corresponding  decrease  in  the  first  mode  mean 
equare  (as  the  vertical  mass  moves  upward) . 
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Fig.  (9a)  Time  history  response  of  second  mode 
excitation  under  excitation  spectral 
density  of  0.022  g2/Hz. 


Fig.  (9b)  Magnification  of  time  history  response  of 
second  mode  excitation  showing  attendance 
of  first  normal  mode. 
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Fig.  (10a)  Mean  square  responses  of  generalized 
coordinates  under  second  normal  mode 
excitation. 
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Fig.  (10b)  Mean  square  responses  of  normal 

coordinates  under  second  normal  mode 
excitation. 
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Fig.  (11a)  Relationship  between  nean  square 

responses  of  generalized  coordinates  and 
excitation  spectral  density  for  various 
values  of  Internal  detuning. 


Fig.  (lib)  Relationship  between  Been  square 

responses  of  normal  coordinates  and 
excitation  spectral  density  for  various 
values  of  internal  detuning. 

(Measurenents  are  taken  for  lower  position  of  the 
upper  aasa . ) 


Figures  (11a)  and  (lib)  show  the  Influence  of  the 
excitation  spectral  density  on  the  normal  mode 
■can  square  responses  of  the  generalized  and  normal 
coordinates,  respectively.  Figure  (lib)  Indicates 
that  the  second  normal  Bode  mean  square  is  relati¬ 
vely  amaller  than  the  flrat  normal  mode  mean  square 
response.  This  suppression  effect  is  due  to  the 
nonlinear  normal  mode  interaction.  However,  the 
saturation  phenomenon,  known  in  deterministic  sys¬ 
tems  with  quadratic  nonlinearity,  la  not  pronounced 
in  the  present  results  since  the  excitation  is  a 
random  process  which  contains  several  frequencies 
each  of  which  may  excite  Che  two  modes.  In  deter¬ 
ministic  excitation,  the  external  and  Internal 
detunings  are  very  important  in  establishing  the 
saturation  phenomenon. 


111,2,3  TJ<}  HOPE- BaiAIlSP 

The  purpose  of  these  tests  is  to  explore  the  beha¬ 
vior  of  the  system  under  random  excitation  which 
covers  both  normal  mode  frequencies.  Due  to  the 
shaker  limitation  the  tests  are  conducted  under 
single  excitation  spectral  density  level  -  0.0026 
g^  /Hz.  A  sample  of  the  time  history  response 
record  is  shown  in  fig.  (12)  which  reveals  the 
presence  of  the  two  modes.  The  amplitude  of  oscil¬ 
lation  of  each  beam  depends  on  the  vertical  mass 
location  which  yields  the  same  internal  resonance 
condition.  Figures  (13a)  and  (13b)  show  the  depen¬ 
dence  of  the  mean  square  response  on  the  internal 
detuning  in  terms  of  generalized  and  normal  coordi¬ 
nates,  respectively.  The  full  points  reveal  linear 
response  characteristics  while  the  empty  points 
show  a  nonlinear  interaction  between  the  two  modes. 


IV.  CONCLUSIONS  AND  DISCUSS I W 

The  results  of  an  experimental  investigation  of 
nonlinear  modal  Interaction  in  a  two -degree- of - 
freedom  structural  model  under  random  excitation 
are  reported.  The  model  equations  of  motion  in¬ 
clude  linear  and  nonlinear  inertia  couplings  of  the 
generalized  coordinates.  The  normal  mode  frequen¬ 
cies  Uj  and  of  the  model  are  adjusted  to  meet 
the  internal  resonance  condition  r  -  2.0.  This 
frequency  ratio  is  found  to  exist  at  two  beam 
length  ratios  /i  -  0.49  and  0.71.  At  these 
locations  the  system  response  characteristics  are 
completely  different  when  the  model  is  excited  by  a 
band  limited  random  excitation.  Three  main  series 
of  tests  are  conducted  to  examine  the  system  res¬ 
ponse  behavior  when  the  first  and  aecond  modes  are 
excited  separately  and  when  both  modes  are  excited 
simultaneously . 

When  the  first  normal  mode  is  externally  excited  it 
is  found  that  the  mean  squares  of  the  two  modes  are 
Increasing  monotonically  with  excitation  spectral 
density.  The  response-excitation  relationship  la 
almost  linear  for  small  excitation  levels.  When 
the  two  beams  are  tuned  to  the  exact  internal 
resonance,  the  response -excitation  relationship 
follows  a  continuous  curve.  For  different  internal 
detuning,  the  response  curves  exhibit  a  disconti¬ 
nuity.  This  feature  is  similar  to  th»  determinis¬ 
tic  characteristics  of  the  same  model. 

When  the  aecond  normal  mode  is  externally  excited, 
the  system  response  is  dominated  by  the  second 
normal  mode  up  to  an  excitation  spectral  density 
level  above  which  the  first  normal  mode  attends  and 


s 


deviation  from  theory  is  attributed  to  the  fact 
that  the  experimental  excitation  it  a  band  Halted 
randon  process,  while  in  theory  it  is  represented 
by  a  wide  band  randon  process .  Another  source  of 
deviation  is  that  the  transformation  into  nornal 
coordinates  is  not  exact  since  it  does  not  take 
into  account  the  effect  of  structural  damping.  To 
eliminate  this  problem,  it  is  convenient  to  adopt 
other  models  whose  generalized  and  nornal  coordi¬ 
nates  are  the  same.  With  new  equipment  and  more 
powerful  shakers  the  first  author  is  currently 
undertaking  an  experimental  research  program 
supported  by  the  NSF. 


Fig.  (12)  Time  history  response  of  two  normal  modes 
excitation. 


Fig.  (13a)  Mean  square  responses  of  generalized 
coordinates  under  band- limited  random 
excitation  of  the  two  normal  modes. 
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Structural  dynamics  with  parameter  uncertainties 


R  A  Ibrahim 

Department  of  Mechanical  Engineering,  Texas  Tech  University,  Lubbock  TX  ?9409 

The  treatment  of  structural  parameters  as  random  variables  has  been  the  subject 
of  structural  dynamicists  and  designers  for  many  years.  Several  problems  have 
been  involved  during  the  last  few  decades  and  resulted  in  new  theorems  and 
interesting  phenomena.  This  paper  reviews  a  number  of  topics  pertaining  to 
structural  dynamics  with  parameter  uncertainties.  These  include  direct  problems 
such  as  random  eigenvalues  and  random  responses  of  discrete  and  continuous 
systems.  The  impact  of  these  problems  on  related  areas  of  interest  such  as 
sensitivity  of  structural  performance  to  parameter  variations,  design  optimiza¬ 
tion,  and  reliability  analysis  is  also  addressed.  The  paper  includes  the  results  of 
experimental  investigations,  the  phenomenon  of  normal  modes  localization,  and 
the  effect  of  mistuning  of  turbomachinery  blades  on  their  flutter  and  forced 


response  characteristics. 


1.  INTRODUCTION 

The  concept  of  uncertainty  plays  an  important  role  in  the 
investigation  of  various  engineering  and  physical  chemistry 
problems.  In  fluid  mechanics,  for  example,  the  inaccuracy  of 
measurements  is  called  “uncertainty”  which  differs  from  the 
concept  of  error  (Kline,  1985).  An  error  in  measurement  is  the 
difference  between  the  true  value  and  the  measured  value.  On 
the  other  hand,  an  uncertainty  is  a  possible  value  that  the  errcr 
might  take  on  in  a  given  measurement.  Because  the  uncertainty 
can  take  on  various  values  over  a  range,  it  is  inherently  random. 
In  control  theory,  the  differential  equations  of  control  systems 
often  involve  uncertain  bounded  state  variables.  The  parameters 
of  transfer  functions  of  certain  models  usually  vary  with  a 
certain  degree  of  uncertainty  (Ashworth,  1982).  Thus  a  prob¬ 
abilistic  transfer  function  can  be  defined  with  uncertain  param¬ 
eters  and  can  lie  anywhere  within  the  ranges  which  are  de¬ 
termined  from  simulation  tests.  The  identification  of  uncertain 
parameters  has  recently  been  examined  by  Skowronski  (1981, 
1984). 

Another  class  of  problems  involving  parameter  uncertainties 
is  the  random  heterogeneity  of  real  media  which  possess  proper¬ 
ties  that  are  described  in  a  probabilistic  sense.  More  specifically, 
these  properties  vary  randomly  with  respect  to  time  and  posi¬ 
tion.  and  thus  constitute  a  random  field.  The  theory  of  wave 
propagation  in  random  media  is  very  complicated  and  involves 
partial  differential  equations  whose  coefficients  are  random 
functions  of  space  and  time.  The  difficulty  of  random  wave 
propagation  problems  stems  from  the  fact  that  the  solution  of  a 
linear  partial  differential  equation  depends  nonlinearly  upon 
the  coefficients  (Chernov,  1960;  Frisch,  1968;  Sobczyk.  1985). 

In  physical  chemistry  the  problem  of  determining  the  vibra¬ 
tional  properties  of  randomly  disordered  crystal  lattices  in¬ 
volves  the  calculations  of  the  frequency  spectrum,  electronic 
energy  levels  of  binary  alloys,  thermodynamic  properties  of 
alloys,  isotropic  mixtures,  and  other  solid  state  phenomena.  Of 
particular  importance  is  the  “normal  localization”  or  “confine¬ 
ment"  phenomenon  which  was  first  reported  by  Anderson 
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(1958).  Anderson  showed  that  the  electron  eigenstates  in  a 
disordered  solid  may  become  localized  and  results  in  a  reduc¬ 
tion  of  metallic  conductivity.  In  structural  dynamics  with 
parameter  uncertainties,  irregularities  may  inhibit  the  propa¬ 
gation  of  vibration  within  the  structure  and  the  vibration  modes 
become  localized.  The  similarities  between  the  propagation  of 
vibration  in  an  elastic  system  and  the  conduction  of  electrons  in 
a  solid  is  discussed  by  Hodges  (1982).  Hodges  and  Woodhouse 
(1983).  and  Pierre  et  al  (1986).  Several  problems  in  physics  and 
physical  chemistry  pertaining  to  crystal  lattice  dynamics  were 
reviewed  by  Elliot  et  al  (1974)  and  recently  documented  in  a 
monograph  by  Bottger  (1983). 

In  structural  dynamics,  uncertainties  arise  from  two  main 
sources  (Prasthofer  and  Beadle,  1975).  The  first  is  a  statistical 
one  and  is  due,  for  example,  to  the  stiffness  or  damping 
fluctuations  caused  by  random  variations  in  material  property  s. 
randomness  in  boundary  conditions,  and  variations  caused  bv 
manufacturing  and  assembly  techniques  The  second  is  nonsta- 
tistical  and  is  due,  for  example,  to  the  inaccuracies  and  assump¬ 
tions  introduced  in  the  mathematical  modeling  of  the  structure. 
In  the  first  class  the  mechanical  properties  of  dynamic  systems 
are  subject  to  a  certain  degree  of  uncertainty  because  the 
physical  properties  of  their  elements  are  not  measured  exactly. 
In  addition,  the  physical  properties  can  experience  variations 
with  the  passage  of  time  as  a  result  of  wear  and  tear  or  just 
inherent  deterioration.  These  properties  should  be  modeled  as 
random  variables  with  a  probability  distribution  representing 
the  distribution  of  the  measured  values.  This  modeling  results 
in  random  eigenvalues,  eigenvectors,  and  random  responses  of 
the  system  in  question.  The  analysis  of  random  eigenvalues  and 
eigenvectors  has  been  a  subject  of  several  studies  by  mathema¬ 
ticians  and  engineers  and  will  be  reviewed  in  section  3. 

Figure  1  shows  five  examples  of  structural  systems  involving 
parameter  and  load  uncertainties.  They  include  “almost"  peri¬ 
odic  structures,  similar  component  subsystems,  mulu-span 
beams,  rocket  fins,  and  turbomachinery  rotors.  The  rocket  fins 
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(a)  Disordered  chain  of  random  mass,  spring,  and  damping 
(Soong  and  Bogdanoff,  1963) 


(b)  Disordered  chain  of 
coupled  pendula 

(Hodges,  1982) 


EI(x),  p(x),  A(x) 


£r  yr 


(c)  Disordered  multi¬ 
span  beam 
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(f)  Nonhomogeneous  column 
with  uncertainties  in 
boundary  conditions  and 
axial  load 


FIG.  I.  Examples  of  disordered  systems 


arc  not  usually  identical  in  their  areas  and  each  fin  has  some 
,  misalignment  with  the  rocket  longitudinal  axis.  For  the  case  of 

turbomachinery  rotors,  there  is  always  some  mass  and  stiffness 
eccentricity  in  the  disks.  Parameter  variations  exist  in  disk 
blades  and  result  in  corresponding  variations  in  the  individual 
natural  frequencies  of  the  blades.  This  problem  is  known  as 
mistuning  (Srinivasan.  1984)  which  may  have  a  significant  effect 
on  the  forced  response  amplitude  of  the  blades  and  also  in  the 
»  value  of  the  flow  speed  at  which  flutter  of  the  blades  occurs. 

Other  examples  include  buned  pipelines,  railroad  trackes.  and 
interconnected  girders.  The  uncertainties  in  these  systems  affect 
to  a  large  extent  their  design  and  operating  performance. 

It  should  be  noted  that  parameter  irregularities  may  cause 
significant  changes  in  the  dynamic  characteristics  of  structural 


systems.  In  particular,  they  may  cause  the  occurrence  of  mode 
localization  which  can  be  used  as  a  means  of  passive  control  of 
vibrations.  In  civil  engineering  the  mechanical  and  strength 
properties  of  the  material  vary  from  one  point  to  another  point 
and  are  seldom  prone  to  certain  in  situ  measurements  but  only 
to  indirect  estimates  (Augusti  et  al,  1984)  The  uncertaintv  of 
these  properties  has  a  direct  relationship  to  the  reliability  of 
such  structures  These  uncertainties  are  usually  manifested  in 
the  applied  loads,  stiffness,  and  theoretical  models  that  are  used 
to  describe  and  relate  loading  and  resistance  The  design  of 
structures  under  conditions  of  uncertainty  implies  a  balancing 
decision  between  risk  of  failure  and  cost  or  weight  (Ang  and 
Tang,  1984;  Frangopol,  1986)  The  nsk  is  an  unavoidable 
consideration  for  structural  optimization  problems  It  has  been 
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customary  in  most  reliability  studies  to  measure  the  risk  by  the 
probability  of  failure  (ie,  the  likelihood  of  occurrence  of  some 
specified  limit  state).  On  the  other  hand,  when  restrictions  and 
constraints  of  the  design  are  imprecisely  described,  the  design 
objective  functions  become  fuzzy  (Zadeh.  1965,  1973;  Brown, 
1980;  Brown  and  Yao.  1983).  Recently,  the  fuzzy  set  theory  has 
been  applied  in  multi-objective  fuzzy  optimization  design  of 
ship  grillage  structures  (Gangwu  and  Suming,  1986). 

The  degree  of  sensitivity  of  structures  to  either  deterministic 
design  changes,  or  stochastic  parameter  variations  is  of  great 
importance  to  the  structural  dynamidst.  In  particular,  it  is 
essential  to  determine  if  small  perturbations  can  result  in  sig¬ 
nificant  changes  of  the  free  or  forced  response  amplitudes.  This 
sensitivity  analysis  is  of  great  concern  to  those  who  are  involved 
in  the  control  of  large  flexible  space  structures  (Meirovitch  et  al, 
1983;  Nurre  et  al.  1984).  These  structures  possess  several  modes 
densely  packed  at  low  frequendes.  When  they  are  descretized, 
model  errors  occur  and  the  free  modes  of  vibration  cannot  be 
determined  accurately.  Thus  when  a  control  system  is  designed 
for  natural  frequencies  whose  values  are  assumed  to  be  exact, 
the  model  errors  and  structural  uncertainties  may  deteriorate 
the  performance  of  the  control  loop,  and  may  even  make  the 
system  unstable.  This  problem  results  in  what  is  known  as 
robustness,  ie,  a  control  system  is  termed  robust  if  it  is  rela¬ 
tively  insensitive  to  model  errors  and  structural  uncertainties. 

This  paper  provides  a  review  of  the  recent  theorems  and 
results  pertaining  to  structural  dynamics  with  parameter  uncer¬ 
tainties.  An  early  account  of  the  subject  was  provided  by  Soong 
and  Cozzarelli  (1976).  Three  main  problems  will  be  addressed. 
These  are: 

1.  Random  eigenvalues, 

2.  Random  response  characteristics,  and 

3  Design  optimization  and  reliability. 


3.  RANDOM  EIGENVALUES 


3.1.  Basic  concept  of  random  eigenvalue 

The  value  of  the  natural  frequency  of  simple  single  degree- 
of- freedom  systems  is  given  bv  the  square  root  of  the  stiffness  to 
mass  ratio.  This  value  is  assumed  by  constant  for  identical 
systems.  However,  experiments  have  shown  that  this  value 
varies  randomly  (Mok  and  Murray,  1965)  because  in  reality  the 
physical  properties  of  the  elements  can  neither  be  measured 
exactly  nor  manufactured  exactly.  Thus,  the  eigenvalues  are 
random  variables  whose  statistical  properties  are  determined  by 
the  random  coefficients  of  the  inertia  and  stiffness  terms  of  the 
equations  of  motion.  Consider  for  example  the  natural  frequency 
of  a  simple  mass-spring  system 

A  —  u~  •  k/m . 


the  variation  of  X  due  to  variations  in  stiffness  k  -  ~k  8  k  and 
mass  m  -  m  +  8m,  may  be  expressed  as  a  Taylor  series 
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where  overbar  quantities  refer  to  mean  values  and  A  -  k  /m 
When  the  variations  8m  and  8k  are  random  variables  the 
natural  frequency  will  be  a  random  variable  The  mean  and 
variance  of  A  can  be  evaluated  as  follows 


Before  reviewing  these  three  problems  the  differences  between 
parametric  random  vibration  and  structural  dynamics  with 
parameter  uncertainties  will  be  discussed  first. 

2.  BETW  EEN  PARAMETRICALLY  EXCITED  AND 

DISORDERED  SYSTEMS 

It  is  very  important  to  distinguish  between  two  types  of 
parameter  variaUons  encountered  in  structural  dynamics  The 
first  type  arises  due  to  random  parametric  excitauon  of  systems 
with  esscnuallv  fixed  properties  while  the  second  class  is  inter¬ 
nal  and  is  associated  with  the  system  when  its  parameters  are 
represented  in  a  probabilistic  sense.  In  the  former  case  the 
system  equations  of  mouon  are  stochastic  differential  equations 
with  random  coefficients  represented  by  random  processes 
(Ibrahim.  1985),  while  in  the  latter  case  the  equauons  of  motion 
are  differential  equations  with  random  parameters  represented 
by  random  variables  (Soong.  1973)  The  methods  of  treaung 
dynamic  systems  under  parametric  random  excitations  are  dif¬ 
ferent  from  those  used  in  solving  differential  equations  with 
random  variable  coefficients  Parametric  random  vibration  is 
basically  a  combination  of  the  theory  of  stochastic  processes, 
stochastic  differential  equations,  and  applied  dvnamics  Systems 
with  parameter  uncertainties  (referred  to  in  the  literature  as 
“disordered  systems"),  on  the  other  hand,  involve  boundary- 
value  problem  and  random  field  theory  (Vanmarcke.  1984)  The 
term  “disorder"  has  been  extensively  used  in  the  literature  to 
distinguish  between  the  case  of  random  perturbation  of  the 
system  parameters  (described  by  a  probabilistic  law)  and  the 
case  when  these  parameters  are  perturbed  in  a  deterministic 
sense 


and 
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The  same  is  applied  when  the  mass  moment  of  inertia  is 
included  in  the  equations  of  motion  Collins  and  Thomson 
(1967)  derived  the  statistical  characteristics  of  principal  mo¬ 
ments  of  inertia  and  principal  axes  directions 

Generally.  the  structural  dynamicist  is  interested  in  de¬ 
termining  the  probabil  ty  that  one  or  more  eigenvalues  lie  in  a 
given  range  or  less  than  a  certain  value  (Boyce.  1968)  However, 
the  probabilistic  description  of  the  eigenvalues  and  the  eigen¬ 
vectors  has  been  examined  for  a  limited  and  simple  class  of 
problems  In  most  cases,  it  is  possible  to  calculate  the  statistical 
functions  (such  as  expectations,  variances,  and  covariance  func¬ 
tions)  of  the  eigenvalues  and  eigenvectors 

The  random  eigenvalue  problem  has  been  examined  for  a 
limited  number  of  linear  discrete  and  continuous  sv  stems  The 
treatment  of  these  systems  is  based  on  the  analvsis  of  random 
matrices  and  random  differential  operators  (Scheldt  and  Purkert 
1983)  The  next  subsecuons  will  review  the  methods  and  main 
results  reported  in  the  literature 


random  variables  c;: 


3.2.  Random  eigenvalues  of  discrete  systems 

The  statistics  of  random  eigenvalues  and  eigenvectors  of 
discrete  systems  may  be  determined  by  using  one  of  three  main 
approaches.  These  are  the  transfer  matrix  method,  the  random 
perturbation  method,  and  the  Monte  Carlo  numerical  simu¬ 
lation  algorithm.  The  transfer  matrix  method  (Kemer  1954, 
1956.  Soong.  1962)  utilizes  a  perturbational  expansion  of  the 
random  eigenvalues  in  terms  of  the  random  perturbations  of  the 
system  parameters  The  perturbauon  method  is  based  on  an 
asymptotic  expansion  and  combines  the  ordinary  perturbation 
and  multivariate  statistical  analysis.  The  multivariate  estab¬ 
lishes  the  probability  distributions  of  random  eigenvalues  in 
terms  of  the  distributions  of  the  matrix  coefficients  in  the 
equations  of  motion  The  Monte  Carlo  method,  on  the  other 
hand,  generates  a  random  sample  of  the  system  random  param¬ 
eters  which  are  used  for  computing  numerically  the  eigenvalues 
and  eigenvectors  for  each  set  of  parameters  in  the  sample 
Monte  Carlo  simulations  are  expensive  since  they  require  a 
large  number  of  numerical  solutions  to  define  the  probability 
level  at  the  tails  of  the  distribution.  This  disadvantage  becomes 
evident  when  one  deals  with  large  or  medium  size  systems 
where  numerical  sinulations  become  unrealistic  on  conventional 
digital  computers.  The  first  two  methods  will  be  outlined  in  the 
next  two  sections. 

3.11.  Transfer  matrix  method 

This  method  was  first  developed  for  disordered  periodic 
lattice  systems  by  Kemer  (1954,  1956)  It  was  adopted  by  Soong 
and  BogdanofT  (1963)  to  examine  the  statistics  of  the  random 
eigenvalues  of  disordered  spring-mass  chain  of  ,V  degrees  of 
freedom  of  the  type  shown  in  Figure  1(a)  Basically  the  method 
is  an  extension  of  the  transfer  matrix  developed  originally  for 
free  vibration  of  deterministic  discrete  systems  (Thomson,  1981) 
The  method  transfers  the  displacement  vector  {X],  of  the  j  th 
mass  into  next  mass  displacement  vector  [  X)t. ic 
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{ X }  -  [  I  -  T]  ...,{X} (4) 

where  I  is  the  unit  matrix  and  (I  *  T]  is  the  transfer  matrix  The 
first  displacement  vector  (X),,  is  related  to  the  last  displace¬ 
ment  vector  (X)v  by  the  relationship 


nii-Tj, 


In  order  to  demonstrate  the  method,  a  periodic  disordered 
chain  with  random  masses  and  constant  equal  spnngs  of  stiff¬ 
ness  K  will  be  considered  Let  the  random  mass  be  defined  by 
the  expression 

m  •  m(  1  <  ).  (6) 

where  m  is  the  mean  value  of  the  mass  and  <  is  a  small 
random  variable  with  zero  mean 

The  transfer  matrix  can  be  written  in  the  form 

[I*T1.-[I*T]  ,*[E](.  (7) 

where  [E]  is  a  perturbational  transfer  matrix  which  results 
from  the  random  perturbations  t 

The  characteristic  equation  can  be  established  from  eq  (5) 
The  roots  of  this  equation  arc  the  system  eigenvalues  ui„  In 
order  to  determine  the  statistical  properties  of  the  eigenvalues  it 
is  necessary  to  express  un  in  terms  of  the  random  variables  t 
It  will  be  assumed  that  the  range  over  which  the  values  of  t  are 
distributed  is  small  and  can  be  explained  in  powers  of  the 


Let  the  random  variables  iJ  be  statistically  independent,  identi¬ 
cally  and  normally  distributed  with  zero  mean.  This  means  that 
the  probability  density  function  of  each  is 

/»(«)“— 7f=«p{ -<:/2o,:}.  (9) 

o,v2n 

where  a,:  is  the  variance  of  the  random  variable  t 

From  the  theory  of  random  processes  (Lamng  and  Battin. 
1956).  it  is  known  that  if  the  random  variables  t,  are  indepen¬ 
dent  and  normally  distributed  the  random  eigenvalues  will  be 
normaly  distributed  with  mean  value  and  variance 

These  two  statistical  parameters  provide  the  ele¬ 
ments  of  the  probability  density  function  of  ui„:.  le 

p(  w„  )  -  — ?=exp{  -  (  u„  -  )'/2o-(j  .  ( 10) 
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Figure  2  shows  p(w„)  and  the  standard  deviation  <j_  for  a 
spnng-mass  chain  of  10  degrees  of  freedom  with  o,  -  0  05  It  is 
seen  that  the  randomness  of  the  masses  results  in  a  considerable 
dispersion  in  the  high  frequency  region  The  standard  deviation 
of  the  random  eigenvalues  increases  with  the  standard  deviation 
of  the  mass  perturbations  e,  according  to  the  formula  (Soong. 
1962) 
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3. 2. 2.  Random  perturbation  method 

The  perturbation  method  for  me  deterministic  eigenvalue 
problem  is  well  documented  (Cole,  1968.  Meirovitch.  1980) 
The  method  has  recently  been  extended  for  random  eigenvalues 
by  Scheldt  and  Purkert  (1983)  The  eigenvalues  of  discrete 
systems  are  usuallv  determined  from  the  conservative  pan  of 
the  system  equations  of  motion  whose  eigenvalue  equation  is 
given  in  the  form 

[  K(  s )  -  X M(  s ) ]  {  x  )  -  ( 0 } .  (12) 

where  K(  s )  and  M(  s  )  are  svmmtinc  stiffness  and  mass  matrices, 
respectively  The  elements  of  these  matnees  are  taken  from  the 
entire  sample  space  S.  ie.  j  c  S  A  and  (x),  are  the  /th 
eigenvalue  and  eigenvector,  respectively  The  random  mainces 
K(s)  and  M (j)  can  be  written  as  the  sum  of  deterministic  and 
random  matrices 

K(.t)-K-k(  O. 

M(  i )  -  M  *  m(  i ) .  (13) 

where  k(.s)  and  mn  )  represent  random  fluctuations  in  the 
stiffness  and  mass  matnees.  respectively,  with  zero  means  such 
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FIG.  2.  Probability  density  functions  and  standard  deviation  of  u. 


Alternatively,  the  problem  can  be  stated  by  transforming  eq. 
(12)  into  the  standard  form 


[A-XI]{x}  -  {0}.  (15) 

where  A  is  the  system  dynamic  matrix  which  is  symmetric 
positive  and  has  the  random  perturbational  form 

A(s)-A  +  «(r).  (16) 

The  deterministic  matrix  A  has  the  simple  eigenvalues 

X,  <X;<  •••  <  X„  (17) 

while  the  random  matrix  A<  s )  has  the  random  eigenvalues 

X,(s)<X:(s)<  •••  <  X;(  s ).  (18) 

it  is  clear  that  the  existence  of  the  first  two  moments  of  the 
eigenvalues  X,(r)  is  implied  by  the  existence  of  the  first  two 
moments  of  the  elements  of  A(  s ). 

The  eigenvectors  (x};  are  normalized  by  the  relation 

(x, .*,)-!.  (19) 


vectors  are  n 
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where  X(/  —  X,  -  Xy,  {  Z},2  “  { Z21,  Z;; . Z2„  }T,  and  the  ele¬ 

ments  of  {  Z  }/2  are  given  by  the  expression 
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On  the  other  hand,  the  correlation  relations  of  the  eigenval¬ 
ues  and  eigenvectors  up  to  the  (k  +  l)-th  order  in  the  perturba¬ 
tions  a,,  are 


R\(j-k)  ”  £[X,.X4]  -£[a„au] 


where  (x;,x; )  denotes  the  scalar  (or  inner)  product  of  the  same 
vector  x;,  ie  { x  )f{x)(.  Introducing  the  two  expansions 
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where  (x),  -  (0.0.. . .  ,0.  l'.O _ 0)  is  the  normalized  eigenvec¬ 

tor  associated  with  X(.  X)t(r)  and  {x(s))(A  are  the  contnbu- 
tions  due  to  the  perturbed  elements  of  Ms)  From  the  analytical 
dependence  of  X(  and  (x),  on  the  elements  of  atr).  Scheldt  and 
Purkert  (1973)  showed  that  expansions  (20)  and  (21)  converge 
at  least  for  sufficiently  small  values  of  the  elements  of  Ms)  The 
homogeneous  terms  X/4(r)  and  {x(r)),j  up  to  fourth  order  are 
given  by  Scheidt  and  Purkert  (1983).  These  terms  can  then  be 
used  to  determine  the  expectations  and  correlation  relations  of 
the  random  eigenvalues  and  eigenvectors.  If  the  correlation 
between  the  elements  of  a(r)  -  [*,,)  are  only  given,  then  up  to 
first-order  perturbation  the  means  of  the  eigenvalues  and  eigen- 
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The  analysis  is  called  first  order  perturbation  if  first-order 
terms  in  expansions  (20)  and  (21)  are  retained  and  higher-order 
terms  are  excluded.  It  is  second  order  if  terms  up  to  second 
order  are  kept.  However,  second-order  perturbation  is  tedious 
and  involves  multivariate  statistical  analysis.  Most  of  the 
analyses  reported  in  the  literature  deal  with  the  first-order 
perturbation. 

Problems  involving  a  random  symmetric  matrix  with  multi¬ 
ple  eigenvalues  of  the  unperturbed  matrix  have  been  treated  by 
Scheidt  and  Purkert  (1983).  The  analysis  consists  in  the  formu¬ 
lation  of  a  convergence  condition  for  the  perturbation  expan¬ 
sions. 

Collins  (1967)  and  Collins  and  Thomson  (1969)  considered 
first-order  perturbation  and  derived  the  eigenvalue  and  eigen¬ 
vector  statistics  of  a  multi-degree-of-freedom  system  in  terms  of 
the  covariance  matrix  of  the  system  elements.  With  reference  to 
the  eigenvalue  eq.  (12)  they  showed  that  the  variations  in  the 
mass  and  stiffness  matrices  result  in  the  following  first  order 
variations  in  the  eigenvalue  and  eigenvector,  respectively: 

x.  L  Tr(*/~*y)+  I  7— ('"/-'”/)  +  (25) 
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If  the  elements  of  the  mass  and_stiffness  matrices  of  eq.  (12) 
are  random  variables  with  means  and  m,  and  variances  o? 
and  a^r  then  the  expected  eigenvalues  and  eigenvectors  are  X, 
and  x„.  respectively,  and  the  variance  of  the  eigenvalue  is 

"2  ax  dx 
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and  p(k,,  k,)  is  the  joint  probability  density  function  for  K 
and  K:,  and  p(/  is  the  correlation  coefficient  for  fc,  and  k,. 
Expressions  for  cov(kr  m, )  and  cov(  m  .  m, )  follow  the  same 
format  of  relation  (8). 

For  a  simple  chain  of  equal  springs  and  masses  with  uncor¬ 
related  random  masses  or  with  random  uncorrelated  stiffnesses, 
Collins  and  Thomson  showed  that  the  standard  deviation  of  the 
frequency  is  governed  linearly  with  the  standard  deviations  of 
the  masses  and  stiffnesses.  The  results  were  confirmed  by  an 
independent  Monte  Carlo  simulation  and  were  very  close  to 
those  obtained  earlier  by  Soong  and  Bogdanoff  (1963).  How¬ 
ever.  these  linear  relationships  disappear  when  correlation  exists 
in  the  masses  or  stiffnesses  and  the  eigenvalues  are  not  closely 
spaced  Recently  Pierre  (1985)  considered  two  different  discrete 
systems  and  employed  a  first-order  perturbation  to  solve  for  the 
statistics  of  their  eigenvalues.  The  first  system  is  a  mass-spring 
chain  with  random  mass  and  the  second  is  a  chain  of  coupled 
pendula  with  random  lengths.  His  results  were  found  identical 
to  those  obtained  by  Soong  and  Bogdanoff. 


Schiff  and  Bogdanoff  (1972a,  b)  derived  an  estimator  for  the 
standard  deviation  of  a  natural  frequency  in  terms  of  second- 
order  statistical  properties  of  the  system  parameters.  The  de¬ 
rivation  was  based  on  the  mean  square  approximation  devel¬ 
oped  by  Bogdanoff  (1965,  1966). 

It  may  be  noticed  that  the  statistical  properties  of  random 
eigenvalues  are  usually  based  on  the  assumption  of  normal 
distribution  of  the  system  random  parameters.  However,  for 
correlated  non-Gaussian  parameters  the  analysis  can  be  per¬ 
formed  in  terms  of  another  set  of  Gaussian  random  parameters 
which  are  evaluated  by  using  the  Rosenblatt  (1952)  transforma¬ 
tion.  This  transformation  has  extensively  been  used  in  reliabil¬ 
ity  analysis  when  the  performance  function  is  nonlinear.  This 
issue  wil]  be  addressed  in  detail  in  section  5.1. 


33.  Random  eigenvalues  of  continuous  systems 
3.3.1.  Methods  of  analysis 

Continuous  systems  may  involve  uncertainties  from  two 
main  sources.  These  are  (Boyce  and  Goodwin  1964): 

( i )  Uncertainties  in  the  geometry  and  the  material  proper¬ 
ties.  The  random  variation  in  space  dependent  parame¬ 
ters  results  in  variations  of  the  differential  operators 
governing  the  free  vibrations  of  the  structure. 

(«)  Uncertainties  in  the  support  mechanism  of  the  system 
(or  the  boundary  conditions). 

The  uncertainties  of  the  first  class  constitute  a  random  field. 
According  to  Vanmarcke  (1984)  the  behavior  of  disordered 
systems  is  governed  by  two  general  laws.  The  first  is  a  statement 
of  “conservation  of  uncertainty”  as  measured  by  the  product  of 
the  variance  by  the  scale  of  fluctuation  of  the  property  in  the 
random  field.  The  scale  of  fluctuation  is  taken  as  the  area  under 
the  correlation  function.  This  product  remains  invariant  under 
linear  transformation  that  preserves  the  mean.  The  second  law 
states  that  the  degree  of  disorder  of  a  homogeneous  random 
field,  as  measured  by  the  direction-dependent  bandwidth  mea¬ 
sure.  tends  to  increase  when  a  random  field  is  subjected  to  local 
aggregation. 

For  the  two  classes  of  uncertainties  the  random  eigenvalue 
has  been  determined  for  a  limited  class  of  dynamical  systems. 
These  include  elastic  strings  and  bars  (Boyce.  1962:  Goodwin 
and  Boyce,  1964)  and  elastic  beams  (Boyce  and  Goodwin  1964; 
Bliven  and  Soong.  1969;  Hoshiya  and  Shah,  1971;  Shinozuka 
and  Astill,  1972:  Vaicailis  1974).  Boyce  (1968)  outlined  a  num¬ 
ber  of  techniques  for  determining  the  statistics  of  the  eigenval¬ 
ues  of  systems  described  by  partial  differential  equations  and 
boundary  conditions  involving  uncertainty  in  their  parameters. 
These  differential  equations  are  of  order  2 n  and  usually  written 
in  the  form 

yw(x)-Jtw(x).  (29) 

subject  to  the  boundary  conditions 

*,(w)-0,  1-1.2 . 2  n.  (30) 

where  & ,  Jf .  and  are  differential  operators  (with  respect  to 
the  spatial  coordinate  t )  whose  coefficients  are  random  vari¬ 
ables.  «  (  x )  is  the  displacement  of  the  system  at  x  Equation 
(29)  involves  values  of  h  and  its  first  In  -  1  derivatives  at  the 
end  points  of  the  interval  in  which  solutions  are  sought  The 
eigenvalue  problem  defined  by  eqs.  (29)  and  (30)  is  assumed  to 
be  self-adjoint  and  positive  definite  The  investigation  of  ran¬ 
dom  eigenvalues  has  been  carried  out  via  analytical  or  numeri¬ 
cal  approaches.  The  numerical  methods  include  the  Monte 
Carlo  simulation  and  stochastic  finite  element  methods  The 


analytical  treatment  of  the  random  eigenvalue  problem  of  sys¬ 
tems  described  by  eqs.  (29)  and  (30)  is  outlined  by  Boyce  (1968) 
and  Scheidt  and  Purkert  (1983).  The  mathematical  methods 
which  have  been  used  to  determine  the  statistical  moments  of 
eigenvalues  are  classified  according  to  whether  the  statistical  or 
nonstatistical  part  of  the  analysis  is  performed  first.  One  class 
consists  of  first  expressing  the  solution  in  terms  of  the  system 
parameters,  without  regard  to  whether  these  parameters  are 
random  or  deterministic.  Having  obtained  such  a  solution,  the 
statistical  properties  are  then  determined.  According  to  Keller 
(1962,  1964)  this  approach  is  referred  to  as  “honest”  and  the 
solution  can  be  determined  by  using  one  of  the  following 
techniques  (Boyce,  1968;  Scheidt  and  Purkert,  1983): 

(i)  Perturbation  methods. 

(//)  Variational  methods. 

( iii  >  Asymptotic  estimate  methods. 

(iv)  Integral  equation  methods. 

The  “honest”  approach  does  not  provide  an  exact  solution  and 
•the  above  four  methods  are  not  suitable  for  every  problem.  For 
example,  the  variational  methods  are  not  suitable  for  structures 
with  random  boundary  conditions.  Variational  methods  and 
integral  equation  methods  are  limited  because  they  only  lead  to 
statements  for  the  first  eigenvalue  of  the  system.  Moreover,  in 
order  to  apply  the  integral  equation  methods,  very  strong 
conditions  for  the  calculation  of  the  mean  of  the  eigenvalues  are 
required.  Under  certain  conditions  pertaining  to  the  spatial 
correlation  function,  the  asymptotic  methods  and  perturbation 
techniques  lead  to  the  same  results.  The  perturbation  methods 
have  less  restrictions  and  are  extensively  used  in  the  literature. 

The  approach,  on  the  other  hand,  is  called  “dishonest” 
(Boyce,  1967)  if  the  statistics  of  the  eigenvalue  problem  are 
directly  determined  by  performing  averaging  analysis  to  the 
system’s  partial  differential  equation  and  its  associated  boundary 
conditions.  The  statistics  can  be  evaluated  by  using  one  of  the 
following  methods: 

( i )  Iteration  methods. 

( // )  Hierarchy  methods  (Haines.  1965.  1967;  Adotnian. 

1983). 

The  iteration  methods  are  based  on  some  assumptions  for  the 
correlation  relations  in  order  to  solve  the  averaged  integral 
equations  of  the  random  eigenvalue.  The  hierarchy  methods 
take  into  consideration  further  equations  so  that  all  statistical 
functions  in  question  can  be  calculated. 

In  a  series  of  papers  by  Purkert  and  Scheidt  (1977, 1979a, b). 
a  number  of  theorems  pertaining  to  functionals  of  weakly 
correlated  processes  '•ncountered  in  the  eigenvalue  problems, 
boundary  value  problems,  and  initial  value  problems  were 
established.  They  treated  the  stochastic  eigenvalue  problem  for 
ordinary  differential  equations  with  deterministic  boundary 
conditions.  The  coefficients  of  the  differential  operator  were 
independently  weakly  correlated  processes  of  small  correlation 
spatial  length.  They  showed  that  as  the  correlation  length 
becomes  very  small,  the  eigenvalues  and  eigenvectors  possess 
Gaussian  distributions.  This  result  has  recently  been  confirmed 
by  Boyce  and  Xia  (1983).  When  the  random  terms  are  not  small 
the  perturbation  method  is  no  longer  valid  and  the  second  term 
in  the  Hermite-Chebychev  expansion  (Ibrahim.  1985)  of  the 
distribution  function  will  not  vanish.  This  implies  that  the 
distribution  of  the  eigenvalue  will  not  be  normal.  Boyce  and 
Xia  (1985)  obtained  the  upper  bounds  for  the  mean  of  eigenval¬ 
ues  through  a  variational  characterization  of  the  eigenvalues. 
For  stochastic  boundary  value  problems  Linde  (1969)  and 


Boyce  (1966, 1980)  considered  a  Sturm- Liouville  problem  with 
a  stochastic  nonhomogeneous  term.  In  their  recent  monograph 
Sheidt  and  Purkert  (1983)  analyzed  the  moments  of  the  ei¬ 
genvalues  and  mode  shapes  of  random  matrices  and  random 
ordinary  differential  operators.  The  calculations  of  these  mo¬ 
ments  were  based  on  perturbation  expansions,  and  so  required 
the  random  terms  to  be  appropriately  small.  Day  (1980)  de¬ 
veloped  a  number  of  asymptotic  expansions  for  the  random 
eigenvalues  and  eigenvectors  of  continuous  systems. 

The  concept  of  the  Wiener  field,  which  is  obtained  by 
replacing  the  time  variable  of  the  Wiener  process  by  a  space 
coordinate,  was  adopted  by  Wedig  (1976,  1977)  as  a  basic 
model  for  randomly  distributed  loadings  or  imperfections  of 
continuous  structural  systems.  The  solution  of  such  boundary 
value  problems  may  thus  be  described  by  integral  equations 
defined  on  the  Wiener  field  and  thus  possesses  the  Markov 
properties.  Wedig  showed  that  these  integral  equations  may  be 
interpreted  in  the  mean  square  sense  via  the  boundary  and 
eigenvalue  problems  of  elastic  structures  with  random  distrib¬ 
uted  imperfections  or  loadings. 

3.3.2.  Applications 

The  random  eigenvalue  of  a  column  under  axial  force  F. 
shown  in  Figure  1(f),  is  described  by  the  second  order  partial 
differential  equation 


d 2  1d2w(x.t) 


d1w(x,t)  d'w(x.t) 

+  F - 1 -  +A(x) - — “  0 


dx ' 


dr 


(31) 


and  the  boundary  conditions: 


EI(0) 


d~w(  x,  t) 

J7- 


It  — 0  ^1 


dw(  x.  t) 


dx 


■0, 


w(0.r)  -0: 


d2w(x,t)  <9w(j;,r) 

EI(L) - h—l-L  +  K2- 


(32) 


dx- 


dx 


■0. 


w(  L.  t)  “  0. 


where  w(.x.  t)  is  the  lateral  displacement  at  distance  x  and  lime 
i.  L  is  the  length  of  the  column.  EI(  x )  is  the  flexural  stiffness, 
and  pA(  x)  is  the  column  mass  per  unit  length.  AT,  and  X.  are 
the  stiffnesses  of  the  end  spnngs.  For  simple  supports  X,  «  X\ 
=>0.  and  for  fixed  supports  X,  —  X,  —  x 

The  solution  of  eq.  (31)  may  be  expressed  in  the  form 


w(x,i) -£t;(X)exp(iwr).  (33) 

Introducing  the  following  substitutions 
X-x/L. 

H  X)  “  /[l  +  a(  X)j. 

A(X)-A[  1  +a(X)]. 

P-FL'-/ E"I. 

X  »  p.4Z/o>~/ El. 


where  a{  -V)  and  a(  .V)  are  random  variables,  eq  (31)  and  the 
boundary  conditions  (32)  for  mode  j  become: 

{[1  +  a(  .V )](/"(  X)}”  -  pL  "(  X)  -X[l  *a(  X)]U  X)  «  0. 

(34) 

[1  +  a(0)j l/"(0)  -  (  X,L/EI)(.  ”<0)  -0  6  (0)  -0; 

[1  +a(l)]l/“(l)  +(X:L/EI)L'(l)-0:  L(1)*0. 

(35) 


7 


where  a  prime  denotes  differentiation  with  respect  to  A",  and 
subscript  j,  indicating  the  mode  number  in  expansion  (33).  is 
removed 

Hoshiva  and  Shah  (1971)  employed  the  standard  perturba¬ 
tion  analysis  to  determine  the  expected  value  and  variance  of 
the  eigenvalue  of  the  nth  mode  by  using  a  linearized  perturba¬ 
tion  technique.  They  found  that  the  variance  of  the  n  th  natural 
frequency  is  proportional  to  the  variances  of  the  stiffness  coeffi¬ 
cients  at  the  boundaries  and  the  axial  load.  This  linear  relation¬ 
ship  implies  that  the  principle  of  superposition  can  be  applied 
in  a  modified  form  For  the  buckling  case,  ie,  when  A  -  0,  the 
eigenvalue  problem  is  reduced  to  determine  the  statistics  of  the 
buckling  eigenvalue  (Augusti  et  al,  1981,  1984).  Shinozuka  and 
Astill  (1972)  considered  the  case  when  both  Kt  and  K2  are 
random  variables. 

The  natural  frequencies  of  transverse  vibration  of  elastic 
beams  were  analyzed  by  Boyce  and  Goodwin  (1964).  They 
considered  the  geometry  of  the  cross-section  of  the  beam  and  its 
support  mechanism  as  random  variables.  The  statistics  of  the 
eigenvalues  were  determined  by  using  three  different  tech¬ 
niques.  These  were  the  perturbation  method,  the  method  of 
integral  equations,  and  numerical  solution.  Bliven  and  Soong 
(1969)  determined  the  statistics  of  the  natural  frequencies  of  a 
simply  supported  elastic  beam  with  random  imperfections  in 
the  beam  stiffness.  The  beam  was  modeled  as  a  lumped- parame¬ 
ter  model  and  the  properties  of  the  frequencies  were  derived  by 
using  a  perturbation  method.  The  stiffness  random  variation 
was  represented  by  the  relation 

EI-El/[l  +  a(x)], 

where  El  is  the  mean  value  of  the  beam  stiffness  and  (x)  is  a 
stationary  random  field  process  with  zero  mean  and  autocorre¬ 
lation  function  given  by  the  relation 

£[a(x,),«(x2)]  «<j2exp(-|*i  -x2[/d),  (36) 

where  d  is  a  non-negative  constant  known  as  the  correlation 
distance. 

The  standard  deviation  of  the  natural  frequency  of  the  beam 
was  obtained  in  the  closed  form 

«„(")  -u(n)o/g(n)  •  (37) 

where  w(n)  is  the  nth  mode  natural  frequency  of  the  uniform 
beam  —  n2v:\ El/mL*  , 

g(n)~f  f'sin2(  nirx,)  sin2  (nrrx, ) 

Jo  Jo 

xexp[  —  | jc,  -  x:l/d]dx,dx:, 

and  m  is  the  beam  mass  per  unit  length. 

Bliven  and  Soong  found  that  when  the  stiffness  fluctuation 
has  zero  correlation  distance  d  »  0,  the  natural  frequency  stan¬ 
dard  deviation  vanishes.  The  standard  deviation  was  found  to 
reach  the  value  of  o^fn)  —  0.5u(n)o  when  the  stiffness  vari¬ 
ation  is  perfectly  correlated  ( d  -*  ao ). 

The  random  eigenvalue  of  a  beam-column  supported  at  its 
ends  by  a  rotary  spnngs  was  examined  by  Shinozuka  and  Astill 
(1972).  The  spring  supports  and  axial  applied  force  were  treated 
as  random  variables.  The  distribution  of  material  and  geometric 
properties  were  considered  correlated  homogeneous  random 
functions.  The  distributions  of  these  properties  were  generated 
by  using  a  Monte  Carlo  simulation  for  multivariate  and  mul¬ 
tidimensional  random  processes  developed  originally  by 
Shinozuka  (1971).  The  mean  and  variance  of  the  eigenvalues 
were  determined  by  using  the  perturbation  analysis  and  Monte 
Carlo  simulation.  It  was  found  that  the  application  of  ap¬ 
proximate  methods,  such  as  the  perturbation  technique  based 
on  exact  or  an  assumed  mode  shape,  causes  a  considerably 


greater  error  for  the  buckling  case  than  in  the  vibration  case. 
Furthermore,  the  perturbation  solution  for  the  eigenvalue  vari¬ 
ance  can  be  approximated  reasonably  well  by  using  an  assumed 
mode  shape  in  place  of  the  unperturbed  mode  shape.  Vaicaitis 
(1974)  employed  a  two- variable  perturbation  expansion  proce¬ 
dure  to  determine  the  eigenvalues  and  normal  modes  of  beams 
with  random  and/or  nonumform  characteristics  which  do  not 
deviate  considerably  from  the  beam  mean  properties.  A  Monte 
Carlo  simulation  was  used  to  determine  the  statistical  averages 
of  beam  eigenvalues  and  mode  shapes.  Two  cases  of  random 
fluctuations  of  beam  cross  section  were  considered.  For  one 
particular  case  there  was  significant  deviation  attributed  mainly 
to  the  fact  that  gradual  change  in  the  beam  stiffness  was 
permitted.  In  this  case  the  beam  is  “soft”  at  one  end  and 
“hard”  at  the  other  end. 

Hart  and  Collins  (1970),  Collins  et  al  (1971),  and  Hasselman 
and  Hart  (1971, 1972)  developed  a  numerical  method  for  com¬ 
puting  the  variance  of  structural  dynamic  mode  properties  by 
using  component  mode  synthesis  which  was  formulated  origi¬ 
nally  by  Hurty  (1964,  1965).  Numerical  solution  provided  rea¬ 
sonable  results  for  lower  modes  even  when  a  relatively  small 
percentage  of  available  component  modes  is  used.  Hart  (1973) 
developed  a  general  algorithm  for  calculating  the  statistics  of 
the  natural  frequencies  and  mode  shapes  of  structures  acted 
upon  by  an  external  static  loading.  This  type  of  problems 
involve  considerable  calculations  due  to  the  fact  that  the  pro¬ 
portionate  axial  load  in  each  member  of  the  structure  is  depen¬ 
dent  upon  the  structural  parameters  which  are  random  vari¬ 
ables.  For  the  two-bar  truss  shown  in  Figure  3  Hart  determined 
the  first  natural  frequency’s  mean  and  standard  deviation.  The 
influence  of  the  static  load  on  the  statistics  of  the  first  natural 
frequency  is  shown  in  Figure  4.  It  is  seen  that  the  standard 
deviation  of  the  natural  frequency  increases  with  the  axial  load. 
The  implication  of  this  increase  was  further  demonstrated  in 
Figure  5  by  using  normal  probability  density  function.  The 
observed  flattening  shape  of  the  probability  density  function 
with  increased  compressive  loading  shows  a  marked  decrease  in 
confidence  with  the  magnitude  of  loading. 

The  random  eigenvalue  problem  of  disordered  periodic  beam 
was  considered  by  Lin  and  Yang  (1974)  They  used  a  first-order 
perturbation  procedure  to  dense  expressions  for  the  vanances 
of  natural  frequencies  and  normal  modes  for  different  cases  of 
random  bending  stiffness  and  span  lengths.  The  natural  fre¬ 
quencies  were  found  to  be  more  sensitive  to  span  sanations 
than  to  bending  stiffness  fluctuation.  It  was  shown  that  if  the 
random  variations  in  bending  stiffness  for  different  spans  are 
uncorrelated  then  there  is  no  effect  on  the  statistics  of  the 
eigenvalues.  The  effect  exists  only  when  there  is  a  correlation  in 
the  random  variation  in  the  individual  spans.  For  a  random 
variation  in  the  span  lengths  it  was  shown  that  the  variance  of 
the  natural  frequency  is  inversely  proportional  to  the  number  of 


FIG  3  Two  bar  truss  (Hart,  1973) 
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FIG.  4.  Variation  in  fundamental  natural  frequency  stausucs  with  applied  load  (Hart,  1973). 


spans.  The  random  imperfections  in  spatial  periodicity  also  with  homogeneous  parameters.  The  standard  deviation  of  the 

resulted  in  variability  in  the  normal  modes.  However,  due  to  the  natural  frequency  of  mth  mode  was  expressed  in  terms  of  the 

arbitrary  choice  of  modal  amplitude  the  variance  of  the  normal  (2m)th  Fourier  coefficients  of  the  random  parameters  and  was 

mode  was  not  a  unique  function  of  space.  represented  as  a  vector  sum  of  their  standard  deviations.  While 

The  statistics  of  natural  frequencies  of  mistuned  blades  of  a  the  normal  modes  of  a  homogeneous  structure  have  a  shape  of 

circumferentially  closed  packet  of  turbomachinery  were  ex-  harmonic  waves  with  symmetrically  located  nodal  diameters, 

amined  by  Ewins  (1973)  and  Huang  (1982).  When  the  bladed  for  a  structure  with  random  parameters  the  mode  shapes  are 

disk  assembly  is  tuned  and  all  the  blades  are  identical  the  complicated  and  the  nodal  diameters  are  located  unsvmmetn- 

natural  frequencies  and  mode  shapes  are  quite  regular.  Each  cally.  It  was  shown  that  these  modes  have  a  shape  involving  not 

mode  may  be  described  as  having  a  particular  number  of  nodal  only  the  main  harmonic,  but  also  an  infinite  number  of  harmon- 

diameters.  just  as  for  an  unbladed  disk.  However,  when  the  ics.  In  addition,  these  random  normal  modes  are  orthogonal 
blades  are  mistuned  to  a  degree  which  might  well  exist  in  despite  their  complicated  form.  Another  important  feature  was 
service,  the  mode  shapes  and  frequencies  becomes  irregular  In  that  the  phase  angles  of  random  normal  modes  are  not  arbitrary 
this  case  the  natural  frequencies  of  the  individual  blades  can  be  (as  in  the  case  of  a  homogeneous  structure)  but  are  random 
randomly  different  from  one  another.  This  problem  is  belonging  variables  independent  of  the  initial  conditions 
to  systems  with  periodic  random  parameters  and  such  systems  Recently,  the  stochastic  finite  element  method  has  been  used 
are  modeled  by  a  stiff  ring  supported  by  transverse  springs  with  by  Nakagiri  et  al  (1985)  to  determine  the  uncertain  eigenvalue 

randomly  distributed  stiffness  and  mass  parameters  Huang  of  fiber  reinforced  plastic  (FRB)  laminated  plates  These  corn- 

adopted  an  exponential  form  for  the  auto-  and  cross-correlation  posite  materials  usually  exhibit  anisotropy  and  heterogeneits 

function  of  the  random  structural  parameters.  This  form  was  The  elastic  constants  may  fluctuate  around  the  mean  values  due 

originally  assumed  by  Hoshiya  and  Shah  (1971)  The  analysis  of  to  some  slackness  during  the  manufacturing  process  w  hich 

Huang  was  based  on  a  spectral  analysis  method  He  found  that  causes  spatial  distribution  of  the  volume  fraction  In  addition, 

the  mean  of  the  natural  frequency  of  the  structure  with  random  another  parameter  known  as  the  stacking  sequence  is  usuallv 

parameters  is  identical  to  the  natural  frequency  of  the  structure  used  as  a  major  design  parameter  of  the  FRB  laminated  plates 


e.\ 

/• 
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The  stacking  sequence  (Vinson  and  Chou,  1975)  implies  a  group 
of  parameters  such  as  elastic  constants,  layer  number,  fiber 
orientation,  and  layer  thickness.  Nakagiri  et  al  considered  the 
effect  of  the  fluctuation  of  the  overall  stiffness  due  to  uncertain 
variation  of  the  stacking  sequence.  The  uncertain  stacking  se¬ 
quence  was  treated  as  a  set  of  random  variables  for  the  case  of 
simplv-supported  graphite/epoxv  plates.  It  was  found  that  the 
eigenvalue  is  more  sensitive  to  the  standard  deviation  of  the 
fiber  orientation,  and  the  effect  of  the  stacking  sequence  is  more 
pronounced  for  the  rectangular  plate  than  for  the  square  one. 
3.3.3.  Normal  mode  localization 

Periodic  structures  with  slight  variations  in  their  periodicity 
can  exhibit  a  phenomenon  known  as  normal  mode  localization. 
This  phenomenon  takes  place  in  a  manner  that  vibrational 
energy  injected  into  the  structure  by  an  external  source  cannot 
propagate  to  arbitrarily  large  distances,  but  is  instead  substan¬ 
tially  confined  to  a  region  close  to  the  source.  Hodges  (1982) 
called  this  phenomenon  as  “Anderson  localization"  due  to 
Anderson  (1958)  who  discovered  mode  localization  in  solid 
state  physics  in  an  attempt  to  understand  electrical  conduction 
processes  in  disordered  solids.  The  effect  of  irregularities  has  a 
similar  effect  to  damping  in  that  it  limits  the  propagation  of 
vibrations  at  large  distances  from  the  excitation  source.  This 
effect  is  mainly  caused  by  confinement  of  the  energy  close  to 
the  source,  not  by  dissipation  of  the  energy  as  it  propagates  out. 

The  phenomenon  of  mode  localization  can  be  well  under¬ 
stood  by  using  the  coupled  pendula  example  [Fig.  1(b)]  which 
was  adopted  by  Hodges  (1982).  Hodges  provided  an  excellent 
explanation  of  mode  localization:  If  all  pendula  are  identical  so 
that  their  individual  natural  frequencies  are  precisely  equal, 
then  the  normal  modes  of  oscillation  when  these  pendula  are 
coupled  together  extend  throughout  the  system,  the  amplitude 
of  oscillation  of  each  pendulum  varies  sinusoidally  with  its 
position  in  space.  On  the  other  hand,  if  the  natural  frequency  of 
oscillation  varies  from  pendulum  to  pendulum  in  some  kind  of 
random  fashion,  then  in  the  limit  of  zero  coupling,  normal 
modes  consist  of  oscillation  of  individual  pendula  at  frequen¬ 
cies  equal  to  their  natural  frequencies.  For  small  coupling  the 
normal  modes  remain  localized  close  to  individual  pendula  and 
the  normal  mode  frequencies  approximate  the  natural  frequen¬ 
cies  of  the  pendula.  Thus  for  a  particular  mode  one  pendulum  is 
oscillating  close  to  its  natural  frequency  with  a  large  motion.  Its 
nearest  neighbors,  unlike  the  ordered  system,  are  driven  off 
resonance,  and  since  the  coupling  is  weak  they  respond  with 
much  smaller  amplitudes.  These  neighbors  in  turn  drive  pendula 
further  out  and  so  on.  but  at  each  step  the  driving  force  and 


irregularities  in  shrouded  blades  of  jet  engine  rotors  can  result 
in  a  stabilizing  mechanism  which  is  closely  connected  with  the 
phenomenon  of  mode  localization.  In  the  framework  of  local¬ 
ization  theory,  the  stabilizing  mechanism  is  explained  based  on 
the  fact  that  the  original  monochromatic  flutter  wave  is  scattered 
into  waves  of  different  and  more  stable  wavelengths  and  inter¬ 
blades  phase  angles.  While  the  effect  of  misiumng  between 
turbomachinerv  blades  is  favorable  in  flutter  (see  also  Kaza  and 
Kielb.  1982)  it  can  lead  to  an  increase  in  amplitude  on  at  least 
one  blade  in  forced  vibration  situation  as  will  be  shown  in 
section  4.1.2. 

For  periodic  multispan  beams  Miles  (1956)  showed  that  the 
natural  frequencies  are  clustered  in  an  infinite  number  of  groups, 
or  bands,  with  n  frequencies  in  each  band,  where  n  is  the 
number  of  spans.  If  a  torsionial  spring  is  placed  at  the  n  -  1 
intermediate  support  location,  then  the  width  of  the  frequency 
bands  diminishes  as  the  spring  constant  k  increases.  In  the 
limit  as  the  spring  constant  goes  to  infinity,  the  beam  becomes 
clamped  at  the  constraint  locations  and  the  width  of  the 
frequency  bands  is  reduced  to  zero.  Pierre  et  al  (1986)  estab¬ 
lished  an  internal  coupling  parameter  which  is  equivalent  to  the 
inverse  of  the  torsional  spring  constant  1  /k9.  For  ke  =  0  the 
spans  are  fully  coupled.  For  large  values  of  the  spring  constant 
and  irregular  spacing  between  supports,  a  multispan  beam  can 
be  regarded  as  a  disordered  chain  of  weakly  coupled  subsys¬ 
tems.  PierTe  (1985)  and  PierTe  and  Dowell  (1986)  developed  a 
theoretical  analysis  for  the  mode  localization  phenomenon  and 
indicated  that  the  free  modes  of  vibration  are  susceptible  to 
becoming  localized  and  the  natural  frequencies  of  the  multispan 
beam  are  in  bands  of  small  width  if  the  spring  constant  is  large 
They  proposed  a  general  criterion  stating  that  localization  may 
occur  if  the  width  of  the  frequency  band  of  the  ordered  system 
is  of  the  order  of.  or  smaller  than,  the  spread  in  individual 
natural  frequencies  of  the  disordered  component  systems 

Pierre  et  al  (1986)  determined  the  free  modes  of  transverse 
vibration  of  a  disordered  two-span  beam  by  using  a  Rayleigh- 
Ritz  formulation  with  the  constraint  conditions  enforced  by 
means  of  Lagrange  multipliers.  They  developed  a  modified 
perturbation  method  to  analyze  the  localized  modes.  Figure  7 
shows  the  mode  shapes  for  tuned  and  mistuned  beam  for 
torsional  spring  parameter  c=  1000.  where  c  =  llk^/El.  I  is 
the  length  of  the  beam  and  E  and  I  are  the  Young's  modulus 
and  area  moment  of  inertia  of  the  beam,  respectively.  For  a 
mistuned  beam  it  is  seen  that  mode  localization  is  manifested  in 
that  the  peak  deflection  is  much  larger  in  one  span  than  in  the 
other  one. 


response  tend  to  diminish  in  magnitude.  A  typical  mode  shape 
diagram  is  shown  in  Figure  6.  In  terms  of  forced  oscillations, 
mode  localization  implies  localization  of  the  response  in  the 
vicinity  of  the  excitation  point. 

The  effect  of  mode  localization  was  examined  by  Bendiksen 
(1984a,b)  and  Valero  and  Bendiksen  (1985)  who  showed  that 


FIG  6  Schematic  diagram  of  ihe  site  amplitude  « (i  I  (or  a  local¬ 
ized  normal  mode  (Hodges.  1982) 


4.  RANDOM  RESPONSE 


The  response  of  linear  structural  components  with  uncertain 
parameters  can  be  determined  by  using  standard  techniques 
such  as  the  impulse  and  frequency  response  functions  and 
perturbation  methods,  or  numerical  approaches  such  as  stochas¬ 
tic  finite  methods  and  Monte  Carlo  simulation.  The  results 
reported  in  the  literature  will  be  reviewed  in  the  next  two 
subsections. 


4.1.  Standard  techniques 


4.1.1.  Simple  structural  components 

In  an  attempt  to  examine  certain  aspects  of  the  dynamical 
response  of  statistically  defined  systems.  Chenca  and  Bogdanoff 
(1958)  and  Bogdanoff  and  Chenea  (1961)  considered  a  linear 
single  degree-of-freedom  system  with  independent  discrete  dis¬ 
tributions  in  the  mass,  damping,  and  stiffness  coefficients  The 
analysis  of  Bogdanoff  and  Chenea  was  based  on  a  partial 
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differential  equation  for  the  response  joint  density  function 
(Kozin,  1961).  This  equation  is  known  as  the  Liouville  equation 
(Soong,  1973)  and  is  identical  to  the  Fokker- Planck  equation 
with  zero  diffusion  coefficient.  Small  dispersions  in  the  system 
parameters  were  found  to  result  in  a  considerable  dispersion  in 
the  frequency  response.  The  impulse  response  of  the  same 
system  was  determined  by  using  the  perturbation  method  by 
Chen  and  Soroka  (1973).  They  considered  a  linear  system 
described  by  the  differential  equation 


X+2  SuJ+^X-fii). 


(38) 


where  the  natural  frequency  is  considered  random  ui„  -  u„  + 
tu„,  w„  is  a  constant  and  the  perturbation  u„  is  a  random 
variable  with  zero  mean,  r  is  a  small  perturbational  parameter 
and  f(t)  is  an  impulse  excitation.  Chen  and  Soroka  derived  the 
solution  of  equation  (38)  by  using  a  perturbational  technique. 
Figure  8  shows  a  sample  of  the  time  history  response  curves  for 
damping  ratio  f  -  0.05.  It  is  seen  that  both  the  mean  and 


standard  deviation  of  the  response  amplitude  are  nonstationarv 
and  the  standard  deviation  is  90  degrees  out  of  phase  from  the 
mean.  The  amplitude  of  the  response  standard  deviation  in¬ 
creases  with  time,  and  gradually  dampens  out  after  it  reaches  a 
certain  level.  For  systems  with  a  very  high  natural  frequency, 
the  uncertainty  in  the  natural  frequency  was  found  to  have  very 
small  effect  on  the  response  statistics.  However,  the  effect  is 
significant  if  the  natural  frequency  is  low.  As  the  damping 
factor  decreases,  the  dispersion  from  the  mean  became  substan¬ 
tial. 

The  response  of  multi-degree-of-freedom  systems  with  ran¬ 
dom  parameters  was  examined  by  Soong  and  Bogdanoff  (1963. 
1964)  and  Chen  and  Soroka  (1974).  Soong  and  Bogdanoff 
determined  the  statistics  of  the  impulse  admittance  and 
frequency  response  of  a  linear  chain  with  random  masses  dis¬ 
tributed  in  a  small  range.  Chen  and  Soroka  developed  a  method 
which  relates  the  statistics  of  response  parameters  to  the  statis¬ 
tics  of  the  system  eigenvalues  and  eigenvectors.  They  showed 
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FIG.  8.  Mean  and  standard  deviation  of  impulse  function  l0S(t)  for  {  -  005  (Chen  and  Soroka.  1973) 


that  the  response  statistics  of  disordered  systems  are  higher 
than  those  of  purely  deterministic  systems.  The  instantaneous 
transient  response  statistics  of  an  undamped  linear  multi-de- 
gree-of- freedom  system,  with  random  stiffness,  subjected  to 
arbitrary  but  deterministic  forcing  functions  was  investigated 
by  Prasthofer  and  Beadle  (1975).  For  the  case  of  an  impulsive 
excitation,  they  found  that  the  growth  of  the  response  uncer¬ 
tainty  is  exponential.  As  the  standard  deviation  of  the  stiffness 
increases  the  response  mean  square  increases  rapidly  with  time. 
For  a  multi-degree-of-freedom  system  the  response  decay  rate 
decreases  as  the  correlation  coefficient  between  the  stiffness 
elements  increases.  The  influence  of  damping  uncertainty  on  the 
frequency  response  of  a  linear  multi-degree-of-freedom  system 
was  examined  by  Caravani  and  Thomson  (1973)  They  de¬ 
termined  the  mean  and  standard  deviation  of  the  response  by 
using  a  linearization  technique  and  a  Monte  Carlo  simulation. 
They  pointed  out  that  an  accurate  estimate  of  the  damping 
coefficients  for  lightly  damped  systems,  in  the  neighborhood  of 
a  natural  frequency,  is  very  important  in  determining  the  mean 
and  standard  deviation  of  the  system  response. 

The  means  and  variances  of  the  frequency  response  func¬ 
tions  of  a  disordered  periodic  beam  were  studied  by  Y ang  and 
Lin  (1975).  Two  types  of  excitation  were  considered.  These  were 
a  concentrated  force  (or  moment)  and  a  distributed  force  con- 
vected  at  a  constant  velocity.  It  was  shown  that  the  magnitude 
of  the  statistical  average  of  the  frequency  response  function  can 
be  considerably  greater  than  the  value  computed  without  taking 
into  account  the  random  variation  in  the  span  lengths.  In  the 
neighborhood  of  resonance  frequencies  the  standard  deviation 
of  the  frequency  response  function  becomes  quite  large,  indi¬ 
cating  greater  uncertainty  in  such  regions.  In  the  case  of  con- 
vected  loading  the  use  of  a  perfect  periodic  model  cannot 
account  for  the  response  in  certain  vibration  modes  while  these 
modes  can  be  induces  in  a  disordered  periodic  beam. 

4.1.2.  Misluited  blade  d  daks 

It  has  been  indicated  in  section  3.3.3  that  the  mistuning  of 
turbomachinery  bladed  disks  could  have  beneficial  effect  in  the 
case  of  blade  flutter.  However,  the  effect  is  reversed  in  the  case 
of  forced  vibration  (Whitehead,  1966;  Ewins,  1969:  Stange  and 
MacBain,  1983).  It  is  believed  that  Tobias  and  Arnold  (1957) 
have  made  the  first  attempt  to  understand  the  effect  of  blade 
mistuning  on  the  response  of  stationary  waves  (modes  traveling 


opposite  to  the  direction  of  disk  rotation  so  as  to  appear 
stationary  to  a  fixed  observer).  An  interesting  and  important 
structural  phenomenon  resulting  from  mistuning  is  the  splitting 
of  a  bladed  disk's  diametral  modes  of  vibration  (modes  having 

1.2 _ _  n  nodal  diameters)  into  “twin"  or  "dual"  modes.  The 

presence  of  dual  modes  characteristics  in  a  bladed  disk  can 
significantly  affect  either  or  both  of  its  aeroclastic  stability  and 
resonant  response  characteristics.  Whitehead  (1966)  showed 
that  there  is  an  upper  limit  to  the  effect  of  mistuning  and  is 
given  approximately  by  the  factor  of  (1  +  JV)/2.  where  N  is  the 
number  of  blades  in  the  row.  This  upper  limit  was  obtained 
under  the  assumption  that  the  damping  forces  are  substantially 
less  than  the  aerodynamic  coupling  forces.  Jay  and  Bums 
(1984)  conducted  a  series  of  rotating  and  unrotanng  test  to 
identify  mistuning.  damping,  split  factors  for  various  diametral 
patterns  and  dynamic  strains  signatures  from  resonant  tests  of 
a  shrouded  fan  blade/disk.  System  mode  responses  to  various 
distortion  patterns  were  found  to  involve  standing  waves  and 
traveling  waves. 

A  number  of  lumped  mass  models  of  bladed  disk  assemblies 
have  also  been  used  to  study  the  effects  of  various  blade 
mistune  distributions  on  the  maximum  resonant  response  of  the 
blades  (Wagner,  1967;  Dye  and  Henry.  1969;  El-Bavoumy  and 
Srinivasan,  1975;  MacBain  and  Whaley.  1984).  The  nature  of 
the  lumped  parameter  models  used  in  these  studies  is  such  that 
individual  blade  response  was  studied  in  terms  of  single-  or 
two-degree-of- freedom  blade  modes  whose  vibratory  response 
was  altered  by  mechanical  coupling  via  the  disk  portion  of  the 
models.  Hence,  the  basis  or  starting  point  for  these  lumped 
mass  models  was  the  individual  blade  resonant  frequencies.  The 
results  showed  how  much  greater  or  smaller  the  individual 
blade  response  would  be  for  a  set  of  mistuned  blades  compared 
to  the  response  of  a  tuned  set  of  blades.  For  a  given  mistuning 
distribution  and  excitation,  the  response  of  the  mistuned  set  of 
blades  was  found  to  be  many  times  greater  or  smaller  (de¬ 
pending  upon  the  disk  circumferential  location)  than  the  re¬ 
sponse  of  tuned  blades.  Ewins  and  Han  (1984)  conducted  a 
sereis  of  case  studies  to  examine  the  influence  of  various  param¬ 
eters  on  the  resonant  response  levels  of  individual  blades  on  a 
disk  They  found,  for  the  case  of  a  33-bladed  disk,  that  mistun¬ 
ing  always  increases  the  highest  resonant  response  level  from 
that  experienced  by  a  tuned  system  but  while  some  blades  are 
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more  highly  stressed,  others  suffered  a  lower  level  and  the  mean 
value  is  roughly  constant.  It  was  also  concluded  that  the  highest 
response  is  always  experienced  by  a  blade  of  extreme  mistime. 

Analytical  investigations  of  mistiming  fall  into  three  cate¬ 
gories  (Griffin  and  Hoosac,  1984).  deterministic  ( Dye  and 
Henerv,  1969:  Ewms.  1973;  El-Bavoumy  and  Srinivasan.  1975), 
statistical  (Huang.  1982),  and  combined  and  statistical  ap¬ 
proaches  (Sogliero  and  Srinivasan,  1980:  Kazan  and  Kielb, 
1982;  Muszvnska  et  al.  1981).  Basu  and  Griffin  (1986)  used  a 
deterministic /statistical  approach  and  developed  a  model  in¬ 
volving  aerodynamic  and  structural  interaction  for  studying  the 
effect  of  mistiming  on  bladed  disk  vibration.  They  found  that 
the  mistiming  effect  significantly  decreases  as  the  density  of  the 
gas  flowing  through  the  turbine  is  decreased.  On  the  other  hand 
the  effect  was  found  to  increase  linearly  with  the  number  of 
blades  on  the  disk. 

4.2.  Stochastic  Unite  element  methods 

Recent  developments  of  stochastic  finite  element  methods 
have  promoted  the  analysis  of  structural  dynamics  with  uncer¬ 
tain  parameters  These  techniques  could  be  broadly  classified 
into  statistical  and  nonstatistical  (Liu  et  al.  1985b).  The  statisti¬ 
cal  approach  is  based  on  numerical  simulation  via  Monte  Carlo, 
stratified  sampling,  and  Latin  Hvpercube  sampling.  A  compara¬ 
tive  discussion  of  these  techniques  is  provided  by 
Mckav  et  al  (1979)  All  simulation  methods  require  that  the 
joint  probability  distributions  of  the  excitation  and  random 
parameters  be  available  However,  these  distribuuons  are  sel¬ 
dom  to  be  available.  Instead,  one  usually  may  assume  that  the 
input  random  variables  are  mutually  independent  and  Gaus¬ 
sian.  If  these  random  inputs  are  non-Gaussian  distributed,  one 
may  use  the  Rosenblatt  (1952)  transformation  to  transform 
non-Gaussian  correlated  variables  to  Gaussian  uncorrelated 
ones  Nonstatistical  approaches  include  numerical  integration 
(Liu  et  al,  1985a.  1986),  second  moment  analysis  (Cornell  1972) 
and  stochastic  finite  element  methods  (Nakagin  et  al.  1984;  Liu 
et  al.  1985a.b:  Hisada  and  Nakagin,  1982:  Hisada  et  al,  1983). 
A  major  advantage  of  these  methods  is  that  the  multivariate 
distribution  functions  need  not  to  be  known  but  only  the  first 
two  moments.  Recently  several  stochastic  finite  element  ap¬ 
proaches  have  been  developed  by  Vanmarcke  and  Grigoriu 
(1983).  Liu  et  al  (1985a,b).  Dias  and  Nagtegaal  (1985),  and 
Mori  and  Ukai  (1986)  Linear  problems  in  structural  mechanics 
with  uncertain  parameters  have  been  solved  by  second-moment 
analysis  (Contreras.  1980:  Nakagiri  et  al,  1984). 

Astill  et  al  (1972)  examined  the  problem  of  impact  loading 
of  structures  with  random  geometric  and  material  properties. 
Their  approach  is  a  combination  of  finite  element  method  and  a 
Monte  Carlo  simulation.  For  the  case  of  an  axisvmmetric 
concrete  cylinder  they  assumed  spatial  distributions  of  Young’s 
modulus  and  density  for  each  realization  of  the  test  cylinder. 
Each  test  cylinder  was  subjected  to  the  same  axial  impact 
loading.  The  algorithm  gave  a  sample  of  100  maximum  stress 
intensities  from  which  the  sample  mean  and  standard  deviation 
were  computed.  For  a  certain  intermediate  location  of  the  test 
cylinder  it  was  found  that  the  axial  stress  is  always  different 
from  the  corresponding  stress  in  a  uniform  cylinder. 

Vanmarcke  and  Grigoriu  (1983)  developed  a  stochastic  finite 
element  analysis  for  solving  first-  and  second-order  statistics  of 
the  deflection  of  structural  members  whose  properties  vary 
randomly  along  their  axis.  The  covariance  matrix  of  these 
element  averages  was  obtained  by  simple  algebraic  operations 
on  the  variance  function  which  in  tum  depends  primarily  on  the 
scale  fluctuation.  Although  this  approach  was  used  to  determine 
the  free  end  deflection  of  elastic  members,  Vanmarcke  and 


Grigoriu  claimed  that  it  can  be  applied  to  determine  the  re¬ 
sponse  statistics  to  external  dynamic  excitations  even  when  the 
statistical  information  about  spatial  variation  of  material  prop¬ 
erties  is  limited.  Recently  Liu  et  al  (1985a. b,  1986)  developed  a 
number  of  probabilistic  finite  elements  methods  for  nonlinear 
structural  dynamics.  These  methods  are  applicable  for  corre¬ 
lated  and  uncorrelated  discrete  random  variables.  For  elastic- 
plastic  bar  with  end  load,  they  (Liu  et  al,  1985b)  computed  the 
mean  and  variance  of  the  displacement  at  the  free  end  by  using 
the  probabilistic  finite  element  and  Monte  Carlo  simulation. 
The  solutions  of  the  two  methods  compared  very  well,  however, 
the  probabilistic  finite  element  approach  required  much  less 
computer  time  than  the  Monte  Carlo  simulation.  Unfortunately 
these  results  did  not  reflect  the  influence  of  parameter  uncer¬ 
tainties  on  the  random  response. 

The  dynamic  response  of  random  parametered  structures 
under  random  excitation  has  been  examined  in  a  number  of 
studies  by  Paez  and  his  group  (Chang,  1985;  Bennett,  1985; 
Branstetter  and  Paez,  1986).  These  studies  provide  computer 
programs  in  a  finite  element  framework  to  establish  response 
moments  on  a  step-by-step  basis.  These  numerical  algorithms 
evaluate  the  system  response  characteristics  at  an  advance  time 
by  using  the  statistical  information  about  response  structural 
characteristics,  and  excitation  at  a  previous  time.  Branstetter 
and  Paez  (1986)  examined  their  computer  programs  for  several 
damped  single  degree  of  freedom  systems  and  several  un¬ 
damped  two  degree-of-freedom  systems.  The  responses  of  these 
systems  to  white  noise  excitations  were  obtained  for  random 
stiffness  parameters  while  all  other  system  parameters  were 
fixed.  It  was  shown  that  single-degree-of-freedom  systems  dis¬ 
play  greater  response  variance  than  systems  with  deterministic 
stiffness.  The  difference  in  response  variance  is  found  to  be 
small  when  the  structure  initial  conditions  are  zero.  The  dif¬ 
ference  increases  and  assumes  an  oscillatory  character  when  the 
initial  conditions  depart  from  zero.  The  mean  response  is  non¬ 
zero  for  structures  with  nonzero  initial  conditions  and/or  non¬ 
zero  mean  load. 

Bennett  (1985)  considered  uncertainties  in  the  stiffness  and 
damping  of  single-  and  multi-degree-of-freedom  structural  sys¬ 
tems.  The  random  variables  of  the  system  parameters  were 
replaced  by  a  deterministic  component  (equal  to  the  mean  of 
the  original  random  variable)  and  a  random  component  with 
zero  mean  and  with  variance  equal  to  that  of  the  original 
random  variable.  For  a  single-degree-of-freedom  system 
Bennett  found  that  the  value  of  the  peak  response  increases 
monotonically  with  the  standard  deviation  of  the  stiffness.  For 
lightly  damped  systems  which  do  not  have  zero  mean,  the 
effects  of  the  damping  randomness  on  the  response  are  less 
pronounced  than  those  obtained  when  the  stiffness  was  random. 
The  standard  deviation  of  the  response  at  the  time  of  peak 
response  was  found  to  increase  with  the  correlation  between  the 
stiffness  and  damping. 

5.  DESIGN  OPTIMIZATION  AND  RELIABILITY 
5.1.  Reliability-based  design 

The  study  of  response  of  disordered  systems  is  very  im¬ 
portant  for  design  purposes.  These  responses  can  help  the 
designer  to  establish  acceptable  tolerances  on  system  compo¬ 
nents.  The  main  problem  which  concerns  the  designer  is  how  to 
govern  the  fluctuations  of  the  system  parameters  for  safe  oper¬ 
ations.  For  example  when  the  values  of  the  elastic  displacement 
of  a  structure  are  significant,  the  problem  is  to  set  up  an 
optimum  standard  of  manufacturing  the  structure  components. 
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Here  the  permissible  fluctuation  in  the  structure  parameters 
becomes  a  restrictive  condition.  Generally,  design  optimization 
of  structures  subject  to  reliability  requirements  is  regarded  as 
the  ultimate  goal  of  any  design  procedure.  The  basic  approach 
in  most  reliability-based  structural  optimization  is  to  impose  a 
set  of  constraints  on  overall  system  reliability  or  probability  of 
failure  (Ang  and  Tang,  1975,  1984;  Moses,  1973;  Parmi  and 
Cohn,  1978).  Another  approach  suggests  to  minimize  the  total 
cost  or  weight  for  a  specified  allowable  overall  failure  probabil¬ 
ity  (Frangopol,  1984a;  Hilton,  1960;  Moses  and  Kinser,  1967). 

One  of  the  main  objectives  of  the  designer  is  to  establish  an 
acceptable  probability  of  failure.  Several  procedures  for  the 
analysis  of  probability  of  failure  of  structures  have  been  devel¬ 
oped  (Frangopol,  1984a,  b,  1985a,  b;  Frangopol  and  Nakib, 
1986;  Kam,  1986;  Moses,  1974;  Moses  and  Kinser,  1967;  Moses 
and  Stevenson,  1970).  In  order  to  establish  a  probability  of 
failure  consider  a  structural  system  subjected  to  a  number  of 
external  loads.  The  structure  is  said  to  survive  if  the  applied 
stress  aD  in  the  built-in  section  due  to  all  external  loads  is 
smaller  than  an  ultimate  limit  stress  ov 

°aSV  (39) 

The  equality  sign  in  eq.  (39)  corresponds  to  the  state  of  the 
collapse  threshold  of  the  structure.  In  general,  for  each  limit 
state,  it  is  possible  to  establish  a  critical  inequality  similar  to  eq. 
(39)  and  identify,  in  the  space  of  the  relevant  parameters,  a 
“safe  region  Sf  (or  success  region)”,  where  the  critical  in¬ 
equality  holds,  and  unsafe  region  &  (or  failure  region),  where  it 
does  not  hold.  These  regions  are  shown  in  Figure  9(a)  according 
to  Augusti  et  al  (1984),  where 

S-aD  and  R  -  o, .  (40) 

In  most  cases  the  applied  load  S  *  S(t)  is  a  random  process, 
while  the  resistance  R,  which  is  calculated  or  measured,  is  a 
random  variable.  For  each  actual  structure,  the  resistance  takes 
up  a  constant  value  R0,  although  uncertain,  and  the  representa¬ 
tive  point  (R,  S)  moves  in  time  up  and  down  the  solid  line  in 
Figure  9(a).  Figure  9(b)  shows  a  possible  realization  of  the 
random  loading  process  S'(t).  The  limit  state  is  attained  when 
S(r)  violates  the  threshold  R0.  The  time  to  failure  zfa>)  can  be 
used  as  a  measure  of  the  structure  reliability.  Alternatively,  one 
can  consider  a  time  interval  (0,  /)  and  then  check  the  critical 
inequality  in  the  worst  possible  condition.  This  can  be  for¬ 
mulated  in  probabilistic  terms  by  stating  that  the  probability  of 
failure  P fail  and  the  complementary  probability  of  success  (reli¬ 
ability)  r  -  PiUC  coincide  respectively  with  the  probability  that 
the  critical  inequality  is  violated  at  least  once  in  the  interval 
(0,  /).  In  space  random  variables,  the  probability  that  a  point 
Q,  which  represents  the  significant  input  and  system  parame¬ 
ters,  falls  either  in  the  failure  region  &  or  in  the  success  region 


y’.  Symbolically,  these  states  are 

Pw-Prob[ecy],  and  Psuc-Prob[{?cy).  (41) 

Among  the  basic  formulations  of  reliability  calculations  are 
the  level  1  and  level  2  approaches.  In  level  1  one  simply  applies 
the  characteristic  safety  factor  y  -  R/S.  In  level  2  one  needs  to 
determine  a  reliability  index  ft  which  measures,  in  units  of 
standard  deviation,  the  distance  between  the  average  point  and 
the  boundary  of  failure  region.  This  means  that  larger  values  of 
f)  imply  smaller  probability  of  failure.  The  probability  of  failure 
is  found  (August  et  al,  1984)  to  be  less  dependent  on  the 
coefficient  fo  variation  8,  —  o(S)/£[S]  of  external  excitation  if 
the  corresponding  coefficient  of  resistance  SR  -  o(R)/E[R]  is 
relatively  large,  where  6{S)  and  S(R)  are  the  standard  devia¬ 
tions  of  the  applied  stress  and  the  resisting  stress,  respectively. 
Level  2  reliability  methods  include  the  estimation  of  the 
mimimum  distance  /3  which  is  regarded  as  a  safety  measure  of 
the  smallest  distance  of  the  surface  separating  the  safe  and 
unsafe  regions  from  the  origin  in  the  space  of  random  vari¬ 
ables  Q. 

Generally  the  level  of  performance  of  any  structural  system 
depends  on  the  properties  of  the  system.  Thus,  it  is  possible  to 
characterize  a  function  g(Q )  known  as  the  performance  func¬ 
tion  such  that 

g(  Q)  >  0  ”  the  safe  state,  and 

(42) 

g(  Q)  <  0  -  the  failure  state. 

Geometrically  the  limit-state  equation  g(Q)  «  0  is  an  n-dimen- 
sional  surface  that  is  referred  to  as  the  “failure  surface.”  The 
performance  function  could  be  linear  or  nonlinear.  The  evalua¬ 
tion  of  the  exact  probability  of  safety  for  nonlinar  performance 
function  is  generally  involved  and  the  determination  of  the 
required  reliability  index  would  not  be  as  simple  as  in  the  linear 
performance  function  (Ang  and  Tang,  1984).  For  correlated 
non-Gaussian  random  variables,  the  safety  index  may  be 
evaluated  in  terms  of  another  set  of  independent  Gaussian 
variables  through  the  Rosenblatt  transformation  (1952). 
Hohenbichler  and  Rackwitz  (1981)  developed  an  algorithm  to 
determine  the  safety  index  by  using  the  Rosenblatt  transforma¬ 
tion. 

Tanaka  and  Onishi  (1980)  developed  a  method  of  regulating 
the  deviations  of  random  parameters  and  derived  a  restrictive 
conditional  formula  in  terms  of  the  permissible  displacement 
(or  natural  frequency)  fluctuation.  The  method  is  based  on  the 
linear  deviation  analysis  with  partial  differential  analysis  to¬ 
gether  with  sequential  linear  programming  (SLP)  for  a  number 
of  restrictive  conditions.  Tanaka  et  al  (1982)  treated  the  optimi¬ 
zation  problem  of  the  allowable  variance  of  random  parameters 
by  using  a  perturbation  method  and  Monte  Carlo  simulation. 
They  computed  the  deviation  of  the  steady  state  response  of 


FIG.  9.  (a)  Safe  and  unsafe  regions  for  (b)  realization  of  S(M  (Augusti  et  al.  1984) 
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structural  systems  involving  random  parameters  with  the  pur¬ 
pose  of  regulating  the  deviation  of  the  random  parameters  when 
the  deviation  of  the  response  is  specified. 

The  techniques  of  mathematical  programming  have  been 
extensively  used  to  minimum-weight  design  of  deterministic 
structures  subject  to  constraints  on  stresses,  displacements,  dy¬ 
namic  response,  and  stiffness  (Moses  and  Kinser.  1967);  Moses 
and  Stevenson,  1970;  Moses  1973,  1974).  The  stochastic  pro¬ 
gramming  of  dynamically  loaded  structures  was  developed 
originally  by  Chames  and  Cooper  (1959)  and  is  well  docu¬ 
mented  by  Rao  (1979).  The  basic  idea  of  this  method  is  to 
convert  the  probabilistic  problem  into  an  equivalent  determinis¬ 
tic  one  by  minimizing  the  expected  value  of  the  objective 
function  subject  to  certain  constraints.  Davidson  et  al  (1977) 
applied  the  mathematical  programming  techniques  for  optimi¬ 
zation  of  structures  subject  to  reliability  requirements.  Their 
work  resulted  in  a  general  formulation  of  the  minimum-weight 
optimization  for  indeterminate  structures  with  random  parame¬ 
ters.  Jozwiak  (1985,  1986)  applied  the  stochastic  programming 
based  on  expected  values  in  the  problem  of  optimization  of 
dynamically  loaded  structures  with  random  parameters.  The 
mean  values  of  joint  displacements  and  their  derivatives  were 
determined  by  solving  the  equations  of  motion  of  the  structure 
under  the  constraints  of  minimum  weight. 

Other  techniques  such  as  multi-objective  optimization  meth¬ 
ods  (Rao,  1982, 1984;  and  Schy  and  Giesy,  1981)  and  fuzzy  sets 
(Zadeh  1965, 1973;  Brown,  1980;  Brown  et  al,  1983)  have  been 
employed  to  the  design  of  simple  structural  elements  and 
aeroplane  control  systems  involving  uncertain  parameters  and 
stochastic  processes.  The  basic  idea  in  multi-objective  design  is 
to  include  all  important  objectives  in  a  vector  objective  func¬ 
tion.  The  problem  of  optimizing  structural  systems  involving 
dynamic  restrictions,  random  parameters,  stochastic  processes, 
and  multi-objectives  has  been  outlined  by  Rao  (1982).  By 
considering  the  imprecision  of  the  restridtons  such  as  use, 
design,  construction,  one  may  assume  that,  some  of  the  con¬ 
straints  and  goals  for  each  of  the  objective  functions  are  fuzzy 
or  imprecise  in  multi-objective  fuzzy  optimization  design.  If  the 
corresponding  expectation  functions  for  objective  and  admis¬ 
sion  for  constraint  are  introduced  it  is  possible  to  quantify  the 
fuzzy  objectives  and  constraints.  Guangwu  and  Suming  (1986) 
employed  the  concept  of  multi-objective  fuzzy  design  optimiza¬ 
tion  for  ship  grillage  structures. 

5.2.  Design  sensitivity  to  parameter  variations 

5.2.1.  Basic  concept  of  sensitivity  analysis 

The  sensitivity  of  a  structural  system  to  variations  of  its 
parameters  is  one  of  the  basic  aspects  in  the  design  of  struc¬ 
tures.  The  sensitivity  theory  is  a  mathematical  problem  which 
investigates  the  change  in  the  system  behavior  due  to  parameter 
variations.  The  basic  concepts  of  sensitivity  theory  are  well 
documented  in  several  books,  see,  eg,  Frank  (1978).  The  sensi¬ 
tivity  problem  can  be  stated  by  defining  the  actual  system 

parameters  represented  by  the  vector  a  -  { a, ,  a2 . a„)T 

which  differ  from  the  nominal  value  a0  by  a  deviation  Aa. 
These  parameters  are  related  to  a  certain  vector  x  which  char¬ 
acterizes  the  dynamic  behavior  of  the  system.  In  structural 
dynamics  the  vector  x  can  be  taken  as  the  system  response 
vector.  The  mathematical  model  of  the  system  response  can  be 
written  in  terms  of  the  first  order  differential  equations 

{*}-{/(x,«,t,F)},  {x(«„)}  -{xn},  (43) 

where  F  represents  the  input  vector. 

Generally  a  unique  relationship  between  the  parameter  vec¬ 
tor  and  the  response  vector  is  assumed.  However,  this  is  not 


possible  in  real  problems  because  they  cannot  be  identified 
exactly.  It  is  a  common  practice  in  sensitivity  theory  to  define  a 
sensitivity  function  S  which  relates  the  elements  of  the  set  of  the 
parameter  deviations  Aa  to  the  elements  of  the  set  of  the 
parameter-induces  errors  of  the  system  function  Ax  by  the 
linear  relationship 

Ax  =  S(a0)Aa.  (44) 

This  relation  is  a  linear  approximation  of  eq.  (43)  and  is  valid 
only  for  small  parameter  variations,  ie,  ||Aa||  «  ll«oll-  S  is  a 
matrix  function  known  as  the  trajectory  sensitivity  matrix  which 
can  be  established  either  by  a  Taylor  series  expansion  or  by 
partial  differentiation  of  the  state  equation  with  respect  to  the 
system  nominal  parameters. 

When  the  system  is  random,  the  function  S  is  referred  to  as 
stochastic  sensitivity  function.  Szopa  (1984)  developed  equa¬ 
tions  for  stochastic  sensitivity  functions  to  determine  the  in¬ 
fluence  of  changes  in  the  initial  conditions  on  the  response. 
These  functions  were  apc,d  to  a  stochastic  nonlinear  oscillator 
with  a  limit  cycle.  It  was  found  that  the  mean  values  and  the 
variances  of  the  stochastic  sensitivity  functions  converge  to 
zero.  Szopa  (1986)  used  the  sensitivity  theory  to  investigate  the 
influence  of  changes  in  system  parameters  on  solutions  of 
dynamical  systems.  The  statistics  of  the  stochastic  sensitivity 
functions  were  found  to  have  finite  values  when  the  response 
exhibit  chaotic  characteristics. 

5.2.2.  Design  derivatives 

Consider  the  eigenvalue  problem  given  by  eq.  (15).  It  will  be 
assumed  that  the  eigenvalues  of  the  system  matrix  A  are 
distinct.  The  elements  of  A  are  function  of  the  system  parame¬ 
ters  a.  The  sensitivity  of  the  free  vibration  of  the  structure  as 
well  as  the  sensitivity  of  its  relative  stability  with  respect  to  any 
parameter  of  A  can  be  characterized  by  the  sensitivity  of  the 
eigenvalues  X,  with  respect  to  the  parameters. 

The  partial  derivative 

(45) 

is  known  as  the  eigenvalue  sensitivity  or  the  eigenvalue  deriva¬ 
tive. 

The  eigenvector  sensitivity  (or  derivative)  of  the  system 
matrix  is  also  given  by  the  partial  differentiation 

S«;  “  dxyda^.  (46) 

The  eigenvalue  sensitivity  has  been  examined  mathemati¬ 
cally  by  McCallev  (1960),  Mantey  (1968),  and  Reddy  (1969). 
Frank  (1978)  developed  a  number  of  formulae  to  determine  the 
eigenvalue  sensitivity.  The  derivatives  of  the  eigenvalues  and 
eigenvectors  are  very  helpful  in  design  optimization  of  struc¬ 
tures  under  dynamic  response  restrictions.  They  have  been 
extensively  used  in  studying  vibratory  systems  with  symmetric 
mass,  damping,  and  stiffness  properties  (Fox  and  Kapoor.  1968; 
Kiefling,  1970)  and  in  nonself-adjoint  systems  (Rogers.  1970: 
Plaut  and  Huseyin.  1973;  Rudisill,  1974).  For  distributed 
parameter  systems,  design  derivatives  of  eigenvalues  were  first 
encountered  in  optimization  studies  by  Haug  and  Rousselet 
(1980)  and  Reiss  (1986).  Reiss  used  a  relatively  simple  method 
to  determine  explicit  results  for  the  design  derivatives  of  eigen¬ 
values  and  eigenvectors.  He  expressed  self-adjoint  operator 
equations  in  terms  of  integral  form  by  using  Green’s  function 
(Reiss.  1983).  Recently  Kuo  and  Wada  (1986)  developed  the 
nonlinear  sensitivity  coefficients  and  correction  terms,  usually 
eliminated  during  the  linearization  process  in  the  Taylor  expan¬ 
sion.  The  nonlinear  correction  terms  were  found  significant  in 
problems  involving  many  finite  element  analyses  where  the  size 
of  the  eigenmatrix  is  of  order  10E06  and  the  difference  in  the 
eigenvectors  may  be  of  order  0.01. 


Lund  (1979)  developed  a  method  to  calculate  the  sensitivity 
of  critical  speeds  of  a  conservative  rotor  to  changes  in  the 
design  using  a  state  vector- transfer  matrix  formulation.  Fritzen 
and  Nordman  (1983)  have  developed  the  eigenvalue  and  eigen¬ 
vector  derivatives  for  general  vibratory  system  (with  nonsym- 
metric  system  matrices)  and  used  them  in  evaluating  stability 
behavior  due  to  parameter  changes  in  rotor  dynamics.  Palazzolo 
et  al  (1983)  presented  a  generalized  receptance  approach  for 
eigensolution  reanalysis  of  rotor  dynamic  systems.  Their  method 
has  the  advantage  of  accomodating  system  modificaion  of  arbi¬ 
trary  magnitude  and  treats  the  modifications  simultaneously. 
Rajan  et  al  (1986)  developed  the  eigenvalue  derivatives  for  the 
damped  natural  frequencies  of  whirl  of  general  linear  rotor 
systems  modeled  by  finite  element  discretization.  For  under- 
damped  modes,  the  eigenvalue  derivative  is  complex.  The  real 
part  represents  the  damping  sensitivity  coefficient  while  the 
imaginary  part  gives  the  whirl  speed  sensitivity.  Rajan  et  al 
showed  that  the  combination  of  design  parameter  and  whir! 
frequency  sensitivity  coefficients  may  be  used  to  evaluate  the 
damped  critical  speed  sensitivity  coefficients. 

In  reliability-based  design  optimization  it  is  useful  to  ex¬ 
amine  the  results  to  sensitivity  analysis  in  order  to  determine 
the  influence  of  the  statistical  parameters  on  the  optimum 
solutions.  The  essential  objectives  of  sensitivity  analysis  of  any 
system  is  to  establish  a  measure  of  the  way  each  response 
quantity  varies  with  changes  in  the  parameters  that  define  the 
system  (Grierson,  1983).  Recently,  Arora  and  Haug  (1979)  and 
Frangopol  (1985a)  have  developed  a  technique  for  determining 
the  reliability-based  optimum  design  sensitivity  of  redundant 
ductile  structures.  Frangopol  investigated  the  sensitivity  of  an 
optimum  design  to  changes  in  the  statistical  parameters  that 
define  the  loading  and  resistance  strength  of  the  structure. 


6.  EXPERIMENTAL  RESULTS 

The  first  attempt  to  measure  the  statistics  of  structural 
modal  frequencies  is  believed  to  be  made  by  Mok  and  Murray 
(1965).  They  carried  out  a  series  of  free  flexural  vibration  tests 
of  a  bar  with  a  stepped  profile  and  a  maximum  variation  in  the 
cross  section  of  50%.  The  predicted  and  measured  results  were 
found  very  close.  Twenty  years  later,  Paez  et  al  (1985,  1986) 
conducted  a  series  of  experimental  investigations  to  measure  the 
random  variation  of  the  natural  frequency  of  a  cantilever  beam. 
One  end  of  the  beam  was  mounted  on  a  fixture  through  a  screw 
and  two  washers,  and  the  other  end  carries  a  concentrated 
mass.  The  torque  in  the  screw  establishes  a  preload  which 
governs  the  stiffness  of  the  beam  at  the  fixture.  Paez  et  al 
conducted  19  experiments  each  with  different  values  of  base 
torque  and  stiffness.  The  variation  of  the  fundamental  frequency 
with  the  base  stiffness  was  obtained  experimentally  and  numeri¬ 
cally  (by  using  a  finite  element  program).  It  was  shown  that  the 
standard  deviation  of  modal  frequency  increases  with  the  mean 
modal  frequency.  Another  interesting  feature  observed  by  Paez 
et  al  was  that  the  magnitude  of  random  variation  in  modal 
frequency  can  become  greater  than  the  spacing  between  modal 
frequencies  as  the  frequency  order  increases. 

The  phenomenon  of  normal  mode  localization  was  first 
examined  experimentally  by  Hodges  and  Woodhouse  (1983). 
Their  model  was  a  thin  high-tensile  steel  wire  stretched  between 
two  supports.  Seven  small  lead  weights  were  securely  attached 
initially  at  equal  lengths  and  then  were  shifted  slightly  to  give  a 
controlled  amount  of  irregularity.  Under  a  step  function  force 
with  repeatable  amplitude  the  string  motion  was  observed  and 
measurements  were  taken  for  the  energy  transmission  from  end 
to  end  of  the  string.  Levels  of  energy  attenuation  in  the  dis- 


FIG.  10.  Comparison  of  experimental  and  theoretical  natural  fre¬ 
quencies  of  the  first  two  modes  for  c  -  281.8  (Pierre  et  al.  1986i 


ordered  case  were  found  in  some  cases  quite  large  (99%)  with 
only  2.4%  standard  deviation  in  the  mass  positions. 

Pierre  et  al  (1986)  conducted  an  experimental  investigation 
to  verify  the  existence  of  localized  modes  for  two  disordered 
two-span  beams  shown  in  Figure  7.  The  beam  was  pinned  at 
both  ends  while  the  third  support  with  variable  torsional  suff- 
ness  was  located  near  the  mid-span.  This  middle  support  can  be 
moved  to  various  locations.  A  pure  excitation  torque  was  ap¬ 
plied  to  the  specimen  beam  near  its  intermediate  support 
Figure  10  shows  the  comparison  between  theoretical  and  experi¬ 
mental  natural  frequencies  of  the  first  two  modes  versus  mistim¬ 
ing  parameter  6/  -  M/I,  where  /  is  the  length  of  the  beam,  and 
A/  is  the  variation  from  the  middle  of  the  beam.  The  coupling 
parameter  c  -  2 ktl/EI,  where  k,  is  the  stiffness  of  the  tor¬ 
sional  spring,  E  and  I  are  the  Young's  modulus  and  the  area 
moment  of  inertia  of  the  beam,  respectively  The  degree  of 
localization  of  a  mode  is  expressed  by  the  ratio  A -A,/ A, 
which  represents  the  peak  deflection  in  one  span  to  the  peak 
deflection  in  the  other  span,  such  that  the  numerator  of  this 
ratio  corresponds  to  the  span  with  smaller  peak  deflection.  This 
peak  ratio  is  shown  in  Figure  11  for  the  two  modes  for  two 
different  values  of  torsional  spring  constant  c.  The  mode  shapes 
of  tuned  and  mistuned  beams  are  shown  in  Figure  12.  It  was 
reported  that  for  $/—  2%  and  c  —  281.8,  the  first  mode  of  the 
mistuned  beam  is  strongly  localized  in  the  second  span,  whereas 
the  one  of  the  tuned  beam  is  collective,  that  is  the  peak 
deflection  is  the  same  in  both  spans. 

A  comprehensive  experimental  and  theoretical  investigations 
were  conducted  by  Ewins  (1976)  to  determine  the  effects  of 
turbomachinery  blades  mistuning.  His  bladed  disk  testpiece 
model  consists  of  24  blades.  A  provision  for  adjusting  the  tune 
of  each  blade  individually  was  accomplished  by  adding  a  num¬ 
ber  of  washers  to  a  nut  and  bolt  attached  near  the  tip  of  each 
blade.  The  test  piece  was  excited  by  placing  an  electromagnet 
close  to  its  surface  and  passing  an  alternating  current  through 
the  magnet.  The  response  of  the  bladed  disk  was  detected  by  a 
set  of  strain  gages  fixed  near  the  root  of  each  blade.  The  natural 
frequencies  were  then  measured  by  adjusting  the  frequency  of 
the  magnet  so  as  to  produce  a  large  response  in  the  strain  gage 
outputs.  The  shape  of  each  mode  was  determined  by  examina¬ 
tion  of  the  relative  amplitudes  of  all  the  blades.  It  was  observed 
that  there  was  a  distinct,  though  complex,  pattern  linking  the 
basic  (tuned)  mode  shape  with  the  mistuned  mode  shape  and 
the  mistuned  pattern,  particularly  for  the  lower  diametral  modes. 
Jay  and  Burns  (1986)  conducted  a  series  of  rotating  and  non- 
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FIG  11  Companion  of  experimental  and  theoretical  peak  ratio  A 


of  the  first  two  modes  for  (ai  c  -  90  4,  <b)  c  -  281  8  (Pierre  et  al.  1986) 
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FIG  12  Measures  first  mode  shape  for  tuned  (a)  and  mistuned  (b) 
two-span  beam  for  c  ■■  281  8  (Pierre  et  al.  1986) 


rotating  tests  to  identify  individual  blade  frequencies,  mode 
shapes,  mistiming,  damping,  and  split  factors  for  diametral 
patterns  of  the  3.  4,  5,  and  6  diametral  mode  families  The  first 
harmonic  of  the  normalized  axiaJ  velocity  defiat  at  the  proper 
mass  flow  rate  was  used  to  construct  a  gust  perturbation  veloc¬ 
ity.  These  spanwise  gust  perturbation  velocities  multiplied  b\ 
tfle  product  of  the  density  and  the  relative  velocity  squared 
results  in  the  normalized  force  parameter  It  was  found  that  am 
increase  in  the  perturbation  force  parameter  results  in  an  in¬ 
crease  in  the  dynamic  stress  in  the  bladed  disk  In  addiuon  the 
perturbation  parameter  does  account  for  the  interaction  be¬ 
tween  the  wake  and  modal  response  of  the  system  as  they  are 
changed  by  aerodynamic  loading. 


7.  CONCLUSIONS 

Several  problems  in  structural  dynamics  involving  parameter 
uncertainties  have  been  treated  in  the  literature  These  problems 
include  the  random  eigenvalue  of  disordered  systems,  normal 
mode  localization,  random  response,  design  opumtzation.  and 
reliability  The  mathematical  theory  of  the  random  eigenvalue 
has  reached  the  maturity  stage,  however,  this  theory  has  not 
been  fully  implemented  to  treat  real  engineering  problems  It  is 
observed  that  some  progress  has  been  made  towards  the  devel¬ 
opment  of  numerical  algorithms  such  as  stochastic  finite  ele¬ 
ment  methods  and  Monte  Carlo  simulations  to  determine  the 
response  of  structural  elements.  These  developments  have  pro¬ 
moted  the  investigation  of  several  problems  including  mistuned 
turbomachinery  bladed  disks,  reliability-based  design  and  de¬ 
rivatives  of  eigenvalues  in  design  optimization.  Few  attempts 
have  been  made  to  employ  new  approaches  such  as  multi¬ 
objective  optimization  and  fuzzy  sets  in  design  optimization 
problems.  It  is  believed  that  these  techniques  will  have  new 
research  avenues  in  many  design  problems.  Another  area  of 
potential  future  research  is  the  optimum  design  sensitivity  in 
reliability-based  design  under  multilevel  reliability  constraints 
to  evaluate  the  significance  of  various  uncertainties  and  ap¬ 
proximations  on  the  optimum  solutions. 

The  problems  treated  in  the  literature  have  been  restricted 
within  the  framework  of  the  unear  theory.  The  limitations  of 
•he  linear  formulation  need  to  be  defined  to  provide  the  struc¬ 
tural  dynamicist  the  influence  of  nonlinearities  as  a  source  of 
uncertainty.  Future  studies  should  include  the  influence  of 
geometric  and  material  nonlinearities.  Experimental  investiga¬ 
tions  are  also  very  important  to  examine  the  influence  of 
parameter  uncertainties  of  composite  structures  on  their  dy¬ 
namic  performance. 
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